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ABSTRACT 

The  design  and  construction,  and  preliminary  experimental  results, 
of  a  test  vehicle  for  studying  plasma  electron  beam  interaction  are 
described  in  associated  Report  A.  Two  connected  chambers  of  the  test 
vehicle,  one  at  hydrogen  pressures  as  high  as  10 “5  Torr  for  generation  of 
plasma  and  the  other  at  pressures  less  than  10 Torr  for  the  generation 
of  a  high-power  electron  beam  in  a  c divergent  flow  gun,  are  described. 

A  demountable  subassembly  for  generating  the  plasma  and  directly 
coupling  S-band  signals  to  and  from  the  plasma  and  the  electron  beam  by 
means  of  a  cavity  is  described.  The  use  of  the  cavities  to  diagnose  the 
plasma  is  described  in  part,  with  a  more  detailed  description  being  pre¬ 
sented  in  associated  Report  B . 

The  electron  gun  used  to  generate  the  5~kv  unity  microperveance  beam, 
having  .20  in.  diameter  in  semi -confined  flow,  is  described.  This  beam  is 
passed  from  the  low-pressure  gun  chamber  into  the  high-pressure  plasma 
chamber. 

A  beam  scanner  for  studying  the  properties  of  the  electron  beam  in 
the  gas  and  plasma  environment  is  also  described.  Use  of  this  scanner  to 
study  the  ion-induced  beam  oscillation  is  described  in  associated  Report  C. 
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I.  INTRODUCTION 


In  attempts  to  achieve  higher  microwave  powers  at  higher 
frequencies,  new  methods  of  harnessing  the  controlled  energies  of 
electron  streams  are  being  investigated.  One  promising  method 
that  has  begun  to  be  studied  in  a  few  laboratories  is  that  of  using 
the  electrons  in  a  plasma  of  a  partially  ionized  low-pressure  gas 
to  interact  with  a  controlled,  linear  electron  beam. 

Use  of  a  plasma  as  an  interaction  "circuit"  in  linear  beam 
tubes  has  advantages  over  a  conventional  circuit.  Fir st, longitudinal 
waves  are  coupled  very  strongly  to  the  electron  beam  even  at  ex¬ 
tremely  high  frequencies.  This  is  due  both  to  the  intimate  mixture- 
possible  of  the  electron  beam  and  the  plasma,  and  to  the  fundamental 
properties  of  the  plasma  at  microwave  frequencies.  This  strong 
coupling  leads  to  very  high  gain  per  centimeter  along  the  patli  of  tin- 
electron  beam  and  should  also  allow  greater  efficiency  in  harnessing 
the  electron  beam.  Second,  for  a  plasma,  there  is  only  a  set  of 
simple  surfaces  removed  from  the  electron  beam  region  necessary 
to  maintain  the  plasma,  whereas  for  slow-wave  circuits  in  con¬ 
ventional  traveling -wave  tubes  at  high  frequencies,  fine,  delicate 
metal  electrodes  are  needed  near  the  electron  beam.  As  a  result, 
much  less  intricate  machining  is  required  with  a  plasma.  There 
is  also  much  less  danger  of  stray  beam  electrons  damaging  the 
electrode  structure  because  of  the  thermal  effects  of  bombardment. 
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Finally,  the  properties  of  the  plasma  can  be  easily  controlled  externally 
by  controlling  the  plasma -generating  discharge  current,  the  gas  pres¬ 
sure,  and  the  magnetic  field  that  is  applied  along  the  direction  of  the  elec¬ 
tron  beam  flow  to  confine  the  electron  beam  and  the  plasma  electrons. 

To  attain  high-power  amplification  of  microwave  signals  using 
the  interaction  of  an  electron  beam  with  a  plasma,  several  problems 
need  to  be  solved.  First,  the  natural  fluctuations  of  the  plasma  and  thus 
the  fluctuation  of  its  microwave  properties  must  be  minimized.  These 
fluctuations  are  functions  of  the  electrode  structure  used  to  maintain 
the  plasma  (as  well  as  of  the  confining  magnetic  field),  the  gas  com¬ 
position  and  pressure,  and  the  plasma  discharge  current.  Reso¬ 
nant  cavities  measure  the  fluctuations  in  the  density  and  loss  of  the 
plasma  as  well  as  microwave  noise.  Next,  the  fluctuations  of  the  elec¬ 
tron  beam  caused  by  its  coupling  to  that  portion  of  the  plasma  that  is 
generated  by  the  electron  beam  must  be  minimized.  This  portion  of 
the  plasma  also  fluctuates,  mainly  because  of  its  interaction  with  the 
electron  beam.  This  beam-generated  portion  of  the  plasma  can 
either  be  generated  in  a  region  where  no  plasma  would  be  generated 
by  any  other  means,  or  it  can  be  generated  as  excess  plasma  in  a 
region  where  moderate  fractional  ionization  has  already  been  brought 
about  by  other  means.  In  either  case,  fluctuations  of  the  beam  and  its 
beam-generated  plasma  can  be  measured  under  controlled,  pulsed 
beam  conditions  using  a  beam-scanning  mechanism  to  diagnose  the 
beam,  with  no  microwave  signals  applied  and  with  the  plasma  discharge 
either  off  or  on.  Next  there  is  the  problem  of  coupling  the  microwave 


signals  through  the  plasma  and  to  the  beam -plasma  system  at  the  core 
of  the  plasma.  This  direct  coupling  is  difficult,  compared  with  coupling 
to  the  electron  beam  before  it  enters  the  plasma,  but  it  is  necessary  to 
use  fully  the  beam-plasma  amplifier's  potential  advantage  in  removing 
the  electrodes  from  close  proximity  to  the  electron  beam.  The  diffi 
culty  arises  because  the  properties  of  the  dense  plasma  are  so  different 
from  a  vacuum,  or  even  a  vacuum  containing  an  electron  beam.  The 
problem  can  be  approached  by  studies  of  the  fringing  fields  of  a  beam- 
plasma  system  of  finite  diameter  and  of  the  interaction  of  such  a  system 
with  a  low-loss  resonant  cavity.  Finally,  there  is  the  problem  of  under¬ 
standing  and  controlling  the  high  small -signal  gain  per  centimeter  of 
beam  travel  in  a  beam-plasma  system,  having  a  finite  geometry  as 
well  as  the  saturating  beam  and  plasma  nonlinearities.  The  small- 
signal  part  cf  this  problem  can  readily  be  solved  by  simultaneous 
analytical  and  experimental  studies,  but  the  large signal  problem  is 
much  more  difficult,  even  when  using  the  best  balance  of  analytical  and 
experimental  studies. 

In  order  to  help  solve  some  of  these  problems,  a  cornpr ehensive 
beam-plasma  tesl  vehicle  was  designed  and  constructed.  This  test 
vehicle  had  two  chambers  connected  by  a  small  cylindrical  tunnel.  A 
pulsed  electron  gun  of  unity  microperveance  and  capable  of  generating 
an  electron  beam  of  up  to  10  kw  was  placed  in  one  chamber.  The 
electron  beam  passed  through  the  tunnel  into  the  other  chamber,  where 
a  flexible  plasma-generating  structure  and  microwave  cavities  for 
coupling  simultaneously  to  the  beam  and  plasma  were  located  in  an 
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internal  subassembly.  Using  hydrogen  gas  at  pressures  of  the  order 

.  3 

of  1  0  Torr  in  the  plasma  chamber,  the  beam -plasma  amplification 
interaction  could  be  investigated  in  detail  at  a  chosen  frequency  of  3,  000 
Me-/  s. 

This  report  includes  a  detailed  description  of  the  main  vacuum 
system,  the  internal  plasma  subassembly,  the  electron  gun  and 
focusing  structure,  and  the  spent  beam-scanning  mechanism.  It  also 
includes  a  discussion  of  the  design,  experimental  methods  used,  and 
preliminary  results. 

Three  reports  are  appended.  These  report  on  the  plasma 
diagnosis  method  and  measurements  ("Plasma  Diagnostics",  by 
D.  M,  Kerr);  on  experimental  results  and  qualitative  descriptions 
of  ion  oscillations  in  linear  beams  ("Ion-induced  Oscillations  of 
Electron  Beams",  by  A.  S.  Gilmour)  and  on  the  theoretical  analysis 
of  small-signal  gain  in  finite  beam-plasma  systems  ("Investigation 
of  the  Interaction  of  an  Electron  Bean  with  a  Plasma",  by  II.  Marantz). 
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II.  VACUUM  SYSTEM 


This  chapter  describes  a  universal  test  vehicle  for  the  beam-plasma 
interaction.  The  vacuum  system  and  other  associated  system  components 
of  this  test  vehicle  are  shown  in  Figures  1  and  2.  The  system  consists 
of  an  electron  gun  region,  which  is  pumped  by  a  400  l/s  Vaclon  pump 
and  which  can  be  sealed  off  from  the  remainder  of  the  test  vehicle  by 
the  vacuum  ball  valve  to  the  left  of  the  gun  region.  To  the  left  of  the 
ball  valve  is  the  plasma  region,  and  to  the  left  of  this  is  the  beam- scan¬ 
ning  area  for  tests  on  the  spent  electron  beam.  Hydrogen  gas  is  provided 
by  a  palladium  leak  tube.  The  leak  rate  is  controlled  by  adjusting  the 
power  applied  to  a  heating  coil  in  the  tube.  In  addition  there  are  the 
following  system  components; 

1.  a  Philips  cold-cathode-discharge  vacuum  gauge, 

2.  an  auxiliary  8  l/s  Vaclon  pump  for  rapid  evacuation  of 
the  test  region, 

3.  an  electromagnet  used  in  the  beam-plasma  region  to 
generate  a  d-c  magnetic  field  up  to  i200  gauss,  and 

4.  a  portable  unit  for  starting  the  Vaclon  pumps. 

A.  SYSTEM  OPERATION 

Operation  of  the  vacuum  system  was  as  follows;  Upon  assembly 

of  the  system,  the  portable  starting  unit  was  used  to  reduce  the  system 

_  3 

pressure  to  at  least  5  x  10  Torr.  This  consisted  of  a  rough  mechanical 
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FIGURE  1.  Vacuum  System  and  Associated  Components  of  Beam-Plasma  Amplifier. 
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pump,  a  cold  trap,  and  VacSorb  pumps.  The  8  l/'s  and  400  l/s  Vaclon 
pumps  were  started  next,  and  the  portable  starting  unit  was  valved  off 
and  disconnected  from  the  system  by  means  of  a  Veeco  quick  coupler; 
then  the  Vaclon  pumps  reduced  the  over-all  system  pressure  to  a  work¬ 
ing  level.  The  plasma  cathodes  were  then  heated,  and  pumping  contin- 

-9 

ued  until  a  steady-state  system  pressure  of  1  x  10  Torr  was  reached 

-7 

in  the  gun  region  and  1  x  10  Torr  in  the  plasma  region.  Next,  the 

8  l/s  Vaclon  pump  was  shut  off,  and  the  hydrogen  leak  source  was 

_  3 

turned  on  to  raise  the  hydrogen  pressure  as  high  as  1  x  10  Torr  in 
the  plasma  region. 

It  should  be  noted  that  in  the  process  of  making  changes  in  the  inter¬ 
nal  subassemblies,  the  ball-valve  arrangement  (shown  in  Figure  3),  in 
conjunction  with  the  valve  in  the  beam-scanner  region  and  that  in  the 
gun  section,  allowed  either  side  of  the  ball  valve  to  be  pumped  or  be  let 
down  to  air  independently  of  the  other  side. 

B.  GAS  FLOW  AND  IONIZATION  CONFIGURATION 

Most  previous  investigators  of  beam-plasma  interaction  phenomena 

have  used  either  mercury  or  cesium  gas  for  the  plasma.  This  is  because 

-4  -  3 

the  mercury  easily  yields  the  pressures  of  interest  in  the  10  —  10  Torr 

region,  in  room-temperature  equilibrium  with  its  liquid  phase,  while 
cesium  is  easily  ionized  by  contact  with  a  heated  tungsten  surface.  Hydro  - 
gen,  however,  is  easily  controlled  in  pressure  1.  y  either  the  controlled- 
Leak  method  or  by  the  heated -generator  method.  Also  since  hydrogen  is 
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FIGURE  3.  Close-up  View  of  Bali -Valve  Section. 
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the  lightest  gas,  its  ions  cause  the  least  sputtering  of  the  cathodes  used 
in  both  the  plasma  discharge  and  the  electron  gun.  Finally,  hydrogen 
does  not  damage  the  cathode  chemically;  in  fact,  it  improves  emission 
somewhat. 

To  avoid  contamination  of  the  cathodes  (especially  the  cathode 

in  the  electron  gun) ,  cold-cathode,  gas-discharge,  getter-ion  (Vaclon) 

pumps  were  chosen  for  this  experiment.  Since  the  pump  makes  it  pos- 

1 

sible  to  pump  hydrogen  faster  than  air,  this  further  adds  to  the  dif¬ 
ference  in  pressure  between  the  high-pressure  plasma  region  (where 
the  hydrogen  source  is  located),  and  the  low-pressure  gun  region  (where 
the  pump  is  located). 

The  400  l/s  Vaclon  pump  was  capable  of  pumping  hydrogen  at 
a  rate  of  over  1000  1/s  at  a  pressure  of  10  ^  Torr  .  Since  the  balL- 
valve  arrangement  between  the  plasma  region  and  the  gun  region  was 
chosen  to  have  a  conductance  of  approximately  10  l/s  fur  hydrogen, 
the  pumping  speed  of  the  over-all  system  for  hydrogen  was  expected 
to  be  approximately  10  l/s;  and  there  should  have  been  a  ratio  of 
100:1  in  pressure  between  the  two  regions  when  the  hydrogen  gas  was 
being  leaked  in.  The  10  ^  pressure  in  the  region  of  the  beam  gun  was 
desirable  (1)  to  minimize  back-ion  bombardment  and  ( Z)  to  minimize 
space -charge  defocusing  and  ion -fluctuation  effects,  which  would  occur 
in  the  gun  region  at  higher  pressures.  Higher  pressure  could  damage 
the  gun  cathode,  during  the  5-kv  negative  pulses,  cause  beam  scalloping 
and  fluctuation,  and  in  addition  prevent  long-life  operation  of  the  400  1/s 
Vaclon  pump. 
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The  pressure  of  10  -  10  Torr  in  the  plasma  region  was  a 

compromise  to  obtain  electron  plasma  densities  sufficiently  high  at  low 
discharge  voltages  without  simultaneously  causing  excessive  scattering 
of  the  electron  beam  and  an  excess  of  H+  ions.  As  a  result,  collisions 
of  beam  electrons  with  the  neutral  gas,  H+  ions  would  be  produced  by 
the  electron  beam  at  a  rate  given  by^ 

n+(t)  =  (n“)  (v^) 

where 

n  (t)  =  the  ion  density  of  the  beam-generated  plasma  (ion/cm  ), 

n^  =  the  electron  density  of  the  beam  (electrons/cm  ), 

_  ^ 

v ^  =  the  ionization  frequency  of  beam  electrons  (sec  ), 

t  =  time  in  seconds  . 

We  can  obtain  from 

vi  =  uQPeip  =  19  x  109p  (sec-1) 

where 

uq  =  the  d-c  velocity  of  electron  beam  (cm/sec)  , 

P  .  =  the  probability  of  ionizing  collisions  per  centimeter  of 

path  per  millimeter  of  mercury  pressure, 

p  -  pressure  in  millimeters  of  mercury. 

We  define  as  the  time  required  for  the  ions  to  produce  space - 

charge  neutralization  of  the  beam  electrons;  this  is  given  by 

„  _  1  .52  x  10'9 

tSCN  "  vi  ~  p 
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The  0.4-a  electron  beam  is  0.20  in.  in  diameter.  These 

assumptions  give  an  electron  density  for  the  beam  of  approximately 

9  3  -4 

10  /cm  .  If  we  now  assume  an  average  pressure  of  .5  x  10  Torr 

in  the  drift-tube  region  and  a  beam-pulse  duration  of  10  (is  ,  we  have 

n+  =  (n")  (M.t) 

+  ,-9  .  ,  3 

n  =10  lons/cm  ; 
max 

i.  e.  ,  space-charge  neutralization  occurs  at  the  end  of  the  pulse. 

It  has  been  shown  (sec  Report  C  by  A.  S.  Gilmour)  that  serious 

beam  fluctuations  of  about  a  few  megacycles  occur  at  the  space-charge 

neutralization  point  during  a  pulse.  This  effect  can  be  delayed  by  lower 

pressure  and  better  ion  drainage. 

3 

Recent  studies  on  ion  defocusing  in  Brillouin-focused  electron 
beams  showed  that  there  was  a  noticeable  distortion  of  the  beam  shape 
when  the  space  charge  of  the  ion  was  only  about  one-fourth  the  space 
charge  of  the  beam.  With  semiconfined  magnetic  focusing,  defocusing 
of  the  space  charge  of  the  beam  by  ions  is  greatly  reduced.  The  radius 
of  the  beam  in  a  gun  with  semiconfined  flow  remains  greater  than  three - 
fourths  of  the  original  space-charge  radius,  whereas  the  radius  of  a 
Brillouin  beam  will,  theoretically,  b  reduced  to  zero  with  complete 
neutralization.  To  avoid  reduction  of  the  beam  radius,  therefore,  a 
gun  with  semiconfined  flow  was  used. 

Further  steps  were  taken  to  minimize  the  effects  of  icn  bom¬ 
bardment  of  the  cathode.  The  profile  of  the  potential  of  the  system  was 
such  that  ions  created  in  the  plasma  region  would  be  accelerated  to 
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the  plasma  cathodes  by  a  potential  of  only  -25  v  .  The  ions  would  be 
reflected  by  an  effective  positive  potential  of  25  v  in  the  grounded 
ball-valve  region.  This  potential  profile  acted  as  an  ion  trap.  Any  of 
the  few  ions  which  did  diffuse  back  through  the  ball-valve  section,  and 
those  created  in  the  high-pressure  ball-valve  region,  would  be  accelerated 
by  the  negative  5-kv  pulses  applied  to  the  gun  cathode,  but  their  bom¬ 
barding  effects  would  be  minimized  because: 

1.  the  electron  gun  has  a  high  (50:1)  area  convergence 
to  localize  the  deterioration  to  the  center  of  the  cathode; 

2.  a  low  (<.001)  duty  cycle  minimizes  the  time  that 
the  ions  bombard  the  cathode. 

C.  HYDROGEN  LEAK  TUBE 

To  maintain  the  desired  pressure  of  hydrogen  gas  in  the  plasma 
region  and  to  control  the  increase  in  pressure  in  this  region  over  that 
in  the  gun  region,  a  hydrogen-leak  apparatus  was  used  in  conjunction 
with  pumping  units.  A  palladium  leak  tube  with  a  maximum  leak  rale 
of  2  l/s  at  1  x  10  ^  Torr  pressure  constructed  in  this  laboratory  is 
shown,  in  Figure  4.  Another  leak  tube  was  purchased,  whose  rate  was 
100  l/s  at  the  same  pressure.  Such  a  tube  was  the  center  of  this  sys¬ 
tem.  A  mixed -source  gas  c.ould  be  used  to  leak  pure,  hydrogen  into  the 
evacuated  system,  because  only  hydrogen  diffused  through  palladium. 

The  theoretical  explanation  for  the  operation  of  the  tube  is  de¬ 
scribed  by  the  diffusion  equation, 
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i 

k  2  -b/T 
=  d  P  6 

where  D  is  the  rate  of  diffusion  (defined  as  the  volume  of  gas  in  cubic 
centimeters  reduced  to  S.  T.  P.  )  that  diffuses  through  one  square  centi¬ 
meter  of  wall  surface  per  second  per  millimeter  of  wall  thickness; 
p  is  the  gas  pressure  in  millimeters  of  mercury;  T  is  the  temperature 
of  the  metal  in  degrees  Kelvin;  d  is  the  thickness  of  the  metal  wall 
in  millimeters;  and  b  and  k  are  constants  for  the  particular  system. 
Calculations  were  made  using  this  equation  with  known  constants  to 
determine  the  surface  area  of  palladium  to  be  used  in  the  Lube  con¬ 
structed  in  this  laboratory  for  a  desirable  range  of  controlling  tempera¬ 
tures. 

Figure  4  gives  an  indication  of  the  size  and  some  indication  of 
the  function  of  the  particular  elements  of  the  tube.  The  hydrogen  leak 
is  a  stainless  steel  tube  through  which  the  hydrogen-nitrogen  gas  mix¬ 
ture  flows.  The  palladium  diffusion  surface  is  simply  a  small  square 
of  metal  brazed  to  a  Kovar  tube  of  proper  size,  which  is  heated  by  a 
filament  that  surrounds  it,  and  which  is  enclosed  in  the  hydrogen  leak 
tube.  A  thermocouple  was  spot-welded  to  the  palladium  surface.  A 
plot  was  made  of  the  characteristic  of  pressure  versus  temperature. 

In  addition,  a  power  supply  for  the  filament  was  constructed,  which 
included  a  variable-voltage  control,  and  a  meter  to  indicate  the  ther¬ 
mocouple  temperature.  To  avoid  the  possibility  of  a  hydrogen  explo¬ 
sion,  a  gas  mixture  of  low  hydrogen  content  is  desirable.  Commercial 
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forming  gas*  (8  percent  hydrogen,  92  percent  nitrogen)  was  used,  and 
since  the  palladium  tube  allowed  only  the  hydrogen  gas  to  diffuse  through 
its  walls,  the  presence  of  nitrogen  in  the  known  mixture  did  not  affect 
the  purity  of  the  hydrogen  plasma  gas. 

D.  FOCUSING  SOLENOID 

To  avoid  radial  diffusion  of  the  plasma  electrons  and  simul¬ 
taneously  focus  the  electron  beam,  a  solenoid  capable  of  establishing 
d-c  axial  magnetic  fields  up  to  1200  gauss,  and  its  associated  power 
supply  (15  a  maximum)  were  constructed.  Provision  was  made  for 
adjusting  the  solenoid  vertically  and  horizontally.  Measurements  of 
the  magnetic  field  of  the  solenoid  alone  and  in  combination  with  the 
solenoid  pole  pieces  indicated  that  a  decrease  in  the  strength  of  the 
magnetic  field  occured  in  the  region  of  the  ball  valve.  It  was  found 
that  the  stem  of  the  ball  Valve,  which  was  machined  at  Cornell  Uni¬ 
versity,  was  made  of  material  that  was  slightly  magnetic,  and  there¬ 
fore  perturbed  the  magnetic  field  enough  to  perturb  the  electron  beam. 
This  perturbation  did  not  seriously  alter  the  beam  performance.  Also, 
since  the  plasma  itself  was  not  located  near  the  ball  valve,  it  was  not 
perturbed  by  this  decrease  in  magnetic  field.  Therefore,  measurements 

Later  in  the  program  it  was  found  desirable  to  have  forming  gas  flow¬ 
ing  when  the  Vaclon  pumps  had  to  be  let  down  to  air.  The  forming 
gas  was  leaked  directly  into  the  400  l/s  Vaclon  pump  (after  shutoff 
and  during  residual  pumping)  before  the  pump  was  let  down  to  air. 

This  made  it  easier  to  restart  the  pump. 
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on  the  various  plasma  parameters  are  to  be  continued  as  scheduled. 

The  high  magnetic  field  in  the  electron  gun  region  was  caused 
by  the  magnet  of  the  400  l/s  Vac  Ion  pump  and  stray  solenoid  fields 
and  was  controlled  by  properly  shielding  the  pump  fields  and  using  a 
small  external  bucking  coil. 

For  optimum  operation  in  future  work,  the  design  that  is  being 
developed  on  another  program  for  the  configuration  of  the  magnetic 
field  for  best  confined  flow  will  be  used.  The  magnetic  field  in  the 
pole  pieces,  plasma  cathode,  and  plasma  regions  agreed  with  theoret¬ 
ical  predictions. 


* 

Linear  Beam  Contract  No,  AF 30(602) - 257 3  i 


-17- 


III.  THE  INTERNAL  PLASMA  SUBASSEMBLY 


A.  MECHANICAL  -  ELECTRICAL  DESIGN 

This  section  describes  the  design  and  construction  of  the  sub- 
assembly  (see  Figures  5  and  6)  to  create  the  Penning-discharge  plasma 
and  to  couple  the  signal  to  and  from  the  beam-plasma  system.  The 
flexibility  of  this  assembly  is  apparent  from  its  demountable  design. 

The  cathode  separation,  interaction  length,  etc.  could  be  changed  by 
changing  the  spacer  dimensions.  The  subassembly  could  be  assembled 
as  a  separate  unit  and  then  mounted  on  main  support  rods  that  were 
heliarced  to  the  flange  of  the  ball  valve.  Because  the  vacuum  system 
and  the  connecting  ball  valve  were  demountable,  the  entire  test  vehicle 
had  great  flexibility. 

The  Penning  discharge  was  generated  by  opposing  impregnated 
ring  cathodes,  which  gave  a  Penning-discharge  geometry  when  com¬ 
bined  with  the  cylindrical  anode  and  longitudinal  magnetic  field.  The 

dual  hot-cathode  Penning  discharge  yielded  an  S-band  density  plasma 

-4 

at  pressures  as  low  as  10  Torr  .  The  electrons  were  trapped  axially 
by  the  electrostatic  field  and  radially  by  the  magnetic  field.  The 
primary  electrons  therefore  had  to  make  several  axial  transits  before 
being  collected  at  the  anode.  The  same  mechanism  trapped  plasma 
electrons  generated  by  primary  electron-neutral  collisions. 

The  anode  structure  was  made  of  three  sections.  The  cylinder 
in  the  interaction  region  between  the  two  cavities  was  of  copper  in  order 
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FIGURE  6.  Internal  Plasma  Subassembly. 


to  reduce  losses  at  microwave  frequencies  and  to  stabilize  the  plasma 
in  that  region.  The  sections  from  the  cavities  to  the  cathodes  were  made 
of  stainless  steel  wires  arranged  in  a  circle  to  allow  observation  of  the 
plasma  and, at  the  same  time,  stabilizing  the  plasma. 

The  ring  cathodes  used  for  the  generation  of  the  plasma  were  of 
the  impregnated  type  and  are  shown  in  Figure  7.  The  cathode  mounting 
is  shown  in  Figure  8a.  The  cathode  was  operated  at  a  negative  potential 
of  approximately  25  v  and  the  anode  at  ground  potential.  Since  the  ball- 
valve  flange  is  at  ground  potential,  there  is  a  possibility  of  backward 
flow  of  electrons  from  the  back  of  the  heated  cathode,  which  could  cause 
arcing.  This  is  prevented  by  the  relatively  cool  molybdenum  backing 
plate  at  cathode  potential,  as  shown  in  Figure  8a.  The  backing  plate 
has  a  hole  to  allow  the  electron  beam  to  pass  and  the  cathode  is  held 
to  this  backing  plate  by  means  of  three  tantalum  wires  120°  apart. 

These  have  small  areas  of  contact  between  the  cathode  and  support 
plate  and  reduce  heat  loss  by  conduction.  A  hollow  stainless- steel 
lube,  polished  on  the  inside,  is  used  as  a  heat  shield  around  the  cathode. 

A  similar  shield  of  a  larger  diameter  is  mounted  to  the  anode  support 
plate,  as  shown  in  Figures  8b  and  c  and  Figure  6. 

The  heatev  leads  were  brought  to  ceramic  standoffs  to  which 
connections  could  be  made.  The  healer  power  recommended  by  the 
manufacturer  was  150  w/in,^  of  cathode  surface,  or  20  w  for  this 
surface.  A  bell-jar  test  was  run  and  the  brightness  temperature  of 
the  cathode  surface  was  measured  as  a  function  of  healer  power.  The 
results  were  plotted  in  Figure  9.  It  was  found  that  because  of  side 


FIGURE  7.  Ring  Cathode  with  Its  Heat  Shield. 
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radiation  losses  and  support  conduction  losses,  heater  power  of 
80  watts  per  heater  was  required  to  attain  the  rated  brightness  tempera¬ 
ture  of  1 190  ° C . 

The  d-c  connections  to  the  subassembly  were  made  to  an  octal 
feedthrough  in  the  ball -val  ve  flange.  The  connection  to  the  feedthrough 
from  the.  heaters,  cathode^tc.  were  made  by  copper  wires  with  ceramic 
insulators  to  prevent  short  circuits. 

R.  PREDICTED  P)  ASMA -DENSITY  PROFILE 

The  radial  density  profile  of  the  electrons  and  ions  in  the  plasma 
generated  by  the  thermionic  electrons  from  a  ring  cathode  can  he  deter¬ 
mined  by  solving  the  continuity  equation, 

---•  l)V‘’u  :  O  ,  (3.1) 

fit 

where  n  is  the  plasma  particle  (primarily  ion)  density  per  unit  volume; 
G  is  the  net  rate  of  particle  creation  per  unit  volume  per  unit  time;  and 
D  is  the  diffusion  constant  Under  steady-stale  conditions  Dn/dt  -  0, 
we  get 

V2n  +  ~  =  0  (3.2) 

For  the  case  of  a  ring  cathode  of  radius  a  and  a  hole  of  radius  b  in  a 
hollow  cylinder  of  radius  R,  we  can  solve  Equation  (3.  2).  In  a  hollow 
cylinder,  n  is  independent  of  <(>  and  z.;  therefore  Equation  (3.2)  re¬ 
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When  the  boundary  conditions  that  n  and  8n/8r  are  continuous  at 
r  =  a,  b  are  used,  the  constants  can  be  evaluated  to  give 

2R  =  0.6,25  in.  ,  2a  =  0.500  in.,  2b  -  0.300  in., 

and  the  normalized  expressions  in  terms  of  C/D  are 

n  2.  =  0.055  ,  0  £  r  <  b  ; 

o 

n5  =  0.161  +  0.073  In  r  -  0.25r2  ,  b  <  r  <  a  ; 

c 

n  •—  =  -0.029  -  0.126  In  r  ,  a  <  r  <  R 

C  _ 

This  profile  is  plotted  in  Figure  10. 

C.  THE  CAVITIES 

Gridless-gap  cavities  of  the  singly  re-entrant  type  were  de¬ 
signed  and  constructed  to  couple  power  on  and  off  the  system.  Copper 
loops  force-fitted  into  the  ends  of  a  1/4  in.  center  conductor  of  rigid 
coaxial  lines  coupled  power  to  the  cavities.  The  copper  center  con¬ 
ductors  were  tapered  along  a  one -wavelength  distance  near  the  loops 
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to  avoid  standing  waves  on  the  line.  The  vacuum  seals  to  the  cavities 
were  made  by  compressing  teflon  discs  at  two  points  of  contact,  i.e.  , 
break  points  in  the  inner  and  outer  conductors  of  the  rigid  coaxial 
transmission  line.  This  seal  was  designed  to  cause  virtually  no  elec¬ 
trical  discontinuity  in  the  transmission  line.  The  vacuum  jacket  seal 
was  made  to  the  outer  conductor  by  a  Veeco  vacuum  coupler  in  the 
output  flange.  The  cavities  were  slightly  overcoupled  to  the  trans¬ 
mission  lines,  since  the  plasma  tended  to  compensate  by  leading  the 
cavities.  The  cavity,  its  coupling  scheme,  and  the  vacuum  seal  are 
explained  in  detail  in  Report  B  of  this  report  and  a  photograph  of  these 
components  is  shown  in  Figure  11. 
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IV.  ELECTRON  GUN  DESIGN 


An  electron  gun  was  designed,  based  on  the  Sperry  Gyroscope 

Company  SM'L-100  gun  used  previously  at  this  laboratory  in  linear 
4 

beam  studies.  The  internal  dimensions  were  the  same  as  those  of  the 
STL-100,  but  the  new  gun  was  made  entirely  of  stainless  steel  and 
set  screws  were  used  in  the  assembly  for  greater  flexibility.  Stain¬ 
less  steel  was  used  instead  of  cold-rolled  steel  so  that  semiconfined 
magnetic  focusing  could  be  used  in  place  of  Brillouin  focusing.  With 
semiconfined  magnetic  focusing,  defocusing  of  the  spao*  cnarge  of  the 
beam  by  ions  is  greatly  reduced.  The  radius  of. the  beam  in  a  gun  with 
semiconfined  flow  will  remain  greater  than  three-fourths  of  the  original 
space-charge  radius,  whereas  the  radius  of  a  Brillouin  beam  will, 
theoretically,  be  reduced  to  zero  with  complete  neutralization.  To 
avoid  reduction  of  the  beam  radius,  therefore,  a  gun  with  semiconfined 
flow  is  preferable. 

Steps  were  taken  to  minimize  the  effects  of  ion  bombardment  of 
the  cathode.  The  profile  of  the  potential  of  the  system  was  such  that 
ions  created  in  the  plasma  region  would  be  accelerated  to  the  plasma 
cathodes  by  a  potential  of  only  about  -25  v  .  The  ions  would  be  reflected 
by  an  effective  positive  potential  of  25  v  in  the  grounded  ball-valve  sec¬ 
tion  between  the  gun  and  the  plasma  region.  This  potential  was  expected 
to  act  as  an  ion  trap.  Ions  which  did  diffuse  back  through  the  ball-valve 
section  and  those  created  in  the  electron  gun  would  have  their  bombarding 
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effects  minimized  by  the  following  means: 


1.  an  electron  gun  with  a  high  (50:1)  area  convergence 
to  localize  the  deterioration  to  the  center  of  the  cathode; 

2.  a  low  duty  cycle  (<.00  1)  to  minimize  the  time  that 
ions  bombarded  the  cathode; 

3.  use  cf  hydrogen  gas,  since  this  ion  causes  less 
cathode  sputtering  than  other  gases  because  of  its  low  mass; 

4.  differential  pumping  resulting  from  the  constriction 
of  the  ball-valve  section  that  permitted  a  100:1  pressure 
differential  between  the  plasma  region  and  the  electron  gun. 

The  specifications  to  be  met  by  the  gun  design  were:  a  5-kv, 

40-a  beam  at  unity  microperveance  with  a  beam  diameter  of  .20  in. 
and  an  area  convergence  of  50:1  .  The  gun  uses  semiconfined  flow 
with  the  focusing  magnetic  field  four  times  that  required  for  Bril- 
louin  focusing.  The  gun  and  focusing  magnet  were  constructed  at 
Cornell.  The  intermediate  details  of  the  gun  are  shown  in  Figure  12. 
The  layout  of  the  gun  assembly  is  shown  in  Figure  13,  and  a  photograph 
of  it  in  Figure  14.  Tiie  supporting  structure  was  designed  both  to  locate 
and  rigidly  support  the  electron  gun.  The  pole  piece  was  one  of  the 
vacuum  flanges  on  the  ball-valve  section  and  the  anode  was  positioned 
on  the  pole  piece  by  a  stainless-steel  ring.  The  gun  was  supported  at 
the  rear  by  another  stainless-steel  ring  and  rods  attached  to  the  rear 
T-section  flange. 

When  the  strength  of  the  magnetic  focusing  field  was  measured 
along  the  axis  of  the  system,  two  problems  arose.  At  the  position  of 
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FIGURE  li.  Details  of  Electron  Gun. 


1/2  Dio.  pumping  holes 


Oxide  cathode 


the  ball  valve,  there  was  a  sudden  dip  in  the  field,  which  was  found  to 
be  caused  by  a  piece  of  magnetic  material  inadvertantly  used  in  the 
construction  of  the  valve;  in  the  gun  region,  the  field  strength  was 
found  to  be  much  higher  than  that  necessary  for  operation  with  semi- 
confined  flow.  It  was  decided  to  leave  the  ball-valve  section  unchanged 
and  construct  a  new  valve  later  if  necessary,  since  the  risk  involved  in 
trying  to  machine  away  the  magnetic  material  was  too  great.  Prelim¬ 
inary  data  could  then  be  taken  despite  the  small  irregularity  in  the  field. 
The  problem  of  the  excess  magnetic  field  in  the  gun  region  was  solved 
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by  making  a  magnetic  shield  that  fitted  around  the  gun,  and,  by  trial 
and  error,  choosing  the  proper  length  to  give  the  desired  ratio  of  the 
field  strength  in  the  plasma  region  to  the  field  strength  in  the  gun  region. 
Figure  15  shows  the  field  shape  before  using  the  shield,  and  Figure  16 
shows  the  field  shape  in  the  final  configuration. 

Final  construction  and  assembly  of  the  electron  gun  proceeded 
routinely  and  no  unexpected  difficulties  arose. 
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IGURE  15.  Magnetic  Forisi.c  r  .eld  on  Axis  of  Beam- 
lasma  System  without  Magnet. i  Shield,  Solenoid  Current 


FIGURE  lb.  Magnetic  Focusing  on  Axis  of  Beam 

Plasma  System  with  Magnetic  Shield  Solenoid  Current 


V.  BEAM-SCANNING  EQUIPMENT 


Since  it  is  planned  to  make  careful  experimental  studies  of  the 
electron  beam  as  it  leaves  the  plasma  interaction  region,  provisions 
were  made  in  the  test  vehicle  to  allow  positioning  of  a  beam-scanning 
mechanism  directly  behind  the  output  plasma  cathode  by  simply  re¬ 
moving  the  copper  collector  now  in  use.  Removal  of  the  collector  was 
accomplished  by  loosening  the  nuts  on  the  ends  of  the  internal  plasma 
subassembly  support  rods  and  main  support  rods,  removing  the  col¬ 
lector  detaching  a  3  in.  extension  section  of  the  support  rods,  and 
replacing  the  nuts. 

It  is  intended  that  the  electron  beam  be  scanned  as  follows: 

1.  With  no  plasma  discharge  and  no  hydrogen  present, 
to  document  the  initial  state  of  the  beam. 

2.  With  hydrogen  present  but  no  plasma  discharge, to 
study  the  beam  fluctuations  resulting  from  beam-generated 
ions;  with  a  plasma  discharge  present  but  with  no  r-f  exci¬ 
tation,  to  study  the  beam  fluctuations  further. 

3.  With  a  plasma  discharge  present  and  with  r-f  exci¬ 
tation,  to  understand  the  gain  saturation  phenomenon. 

In  addition  to  the  results  directly  applicable  to  this  study  information  will 
be  obtained  to  allow  comparison  between  the  beam- scanning  experiments, 
already  completed  by  Gilmour  (see  Report  G  of  this  report),  on  a  Bril- 
louin-focused  beam  partially  defocused  by  ions  and  our  experiments 
under  conditions  of  semieonfined  flow. 
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The  beam-scanning  mechanism  to  he  used  is  machined(but  not 

4 

assembled.  It  is  similar  to  the  scanner  reported  by  Gilmour  and  a 
summary  follows  for  convenience.  In  addition, Figures  21  and  22  in¬ 
dicate  the  type  of  information  that  can  be  derived  from  such  a  scanner. 

BEAM  SCANNING  MECHANISM 

1.  Control  of  Faraday  Cage 
Motion  inside  Vacuum  System 

A  simplified  sketch  of  the  part  of  the  beam-scanning  mechanism 
inside  the  beam  tester  is  shown  in  Figure  17.  Most  of  the  electron  beam 
will  be  collected  by  the  beam-collecting  plate,  which  is  carbonized  to 
reduce  the  secondary  electron  ratio.  A  small  portion,  however,  will 
be  allowed  to  pass  through  the  0.010  in.  aperture  in  the  center  of  the 
plate  to  a  Faraday  cage.  The  electron  current  from  the  cage  will  go 
out  through  a  coaxial  vacuum  feedthrough  in  the  bottom  control  rod  to  an 
indicating  device  outside  the  tester.  The  collecting  plate  will  be  posi¬ 
tioned  in  the  plane  perpendicular  to  the  axis  of  the  solenoid  by  rotations 
of  the  two  control  rods.  As  shown  in  the  sketch,  the  plate  could  be 
moved  horizontally  by  rotating  the  bottom  rod.  Vertical  motion  could 
be  produced  in  a  similar  manner  by  rotating  the  side  rod.  Teflon  bear¬ 
ings,  which  are  used  to  prevent  binding  between  stainless-steel  parts, 
are  attached  to  the  back  of  the  plate  and  slid  on  the  lever  arms  attached 
to  the  side  and  bottom  rods  during  horizontal  and  vertical  motion.  Al¬ 
though  rotational  motions  are  converted  to  translational  motions  in  posi¬ 
tioning  the  cage,  they  could  be  considered  to  be  linear,  since  the  distances 
moved  are  small  compared  to  the  length  of  the  lever  arms. 
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FIGURE  17.  Simplified  Sketch  of  Beam- 
•  canning  Mechanism  inside  Beam  Tester. 


The  collecting  plate  and  positioning  rods  are  illustrated  in  Fig¬ 
ure  18.  The  disk  shown  behind  the  collecting  plate  and  the  lever  arms 
are  used  to  center  the  entire  mechanism  in  the  drift  tube.  As  shown  in 
Figure  19,  this  disk  is  guided  by  a  teflon  bearing  that  rode  on  a  rail  in 
the  drift  tube  and  is  prevented  from  touching  the  drift  tube  by  two  spring- 
loaded  teflon  blocks.  Axial  motion  of  the  cage  assembly  is  produced  by 
moving  the  control  rods  through  the  Veeco  quick  couplings. 
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Spring  -loaded 
teflon  blocks 


FIGURE  19.  Disk  Which  Centers  Faraday  Cage  Assembly. 

2,  Control-Rod  Positioning  Mechanism 
The  angular  position  of  the  control  rods  is  adjusted  by  means  of 
the  micrometers  shown  in  Figure  20.  Each  micrometer  contacts  a  lever 
arm  at  a  radius  that  is  twice  the  radial  position  of  the  cage  inside  the 
tester;  therefore,  incremental  changes  in  the  cage  position  are  deter¬ 
mined  by  taking  one -half  the  incremental  micrometer  readings.  The 
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FIGURE  20.  Control-Rod  Positioning  Mechanism. 


cage  is  positioned  axially  by  using  a  gear-and-rack  arrangement.  In 

4 

the  first  measurements  taken  on  an  electron  beam  by  Gilmour,  he  posi¬ 
tioned  both  micrometers  manually.  More  than  an  hour  was  required  to 
obtain  a  detailed  cross-sectional  plot  of  the  current  densities  in  the  beam 
at  one  axial  position.  Since  it  was  expected  that  hundreds  of  beam  cross 
sections  would  be  needed  to  determine  accurately  the  scallop  wavelength 
of  the  beam  as  a  function  of  the  solenoid  field  and  as  a  function  of  the 
beam  voltage,  a  motor  drive  was  attached  to  the  horizontal  positioning 
micrometer.  In  addition,  the  automatic  reversing  switches  shown  in 
Figure  20  were  installed  to  force  the  motor  to  sweep  the  Faraday  cage 
back  and  forth  automatically.  A  ten-turn  helipot  was  geared  to  the  motor 
to  provide  a  voltage  proportional  to  the  horizontal  position,  and  this 
voltage  was  applied  to  the  x  input  of  an  x-y  recorder.  The  current 
from  the  Faraday  cage  was  fed  to  the  y  input  of  the  recorder.  The 
vertical  position  in  the  beam  was  indicated  in  the  recorder  by  changing 
the  y  reference  position.  A  typical  beam  cross  section  taken  in  this 
manner  is  shown  in  Figure  21, and  Figure  22  shows  the  variation  of  beam 
radius  with  axial  position. 
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VI.  EXPERIMENTAL  RESULTS 


The  generation  and  preliminary  diagnosis  of  the  plasma  in  the 
internal  plasma  subassembly,  as  well  as  the  generation  and  preliminary 
testing  of  the  electron  beam  from  the  semiconfined  flow  convergent  gun, 
have  been  carried  out.  The  microwave  cavities  were  used  to  diagnose 
the  plasma  and  to  test  the  electron  beam. 

The  Penning-discharge  plasma  operated  at  a  nearly  constant 
value  of  40  volts  and  could  be  varied  in  density  by  varying  the  discharge 
current  from  a  few  milliamps  to  a  few  hundred  rnilliamps.  The  plasma 
density  also  could  be  varied  by  controlling  the  hydrogen  pressure,  which 
in  turn  was  controlledby  a  variable  hydrogen  leak,  and  by  the  magnetic 
field.  The  hydrogen  pressure  was  over  100  times  as  great  in  the  plasma 
chamber  as  it  was  in  the  gun  region. 

The  electron  density  and  collision  frequency  are  found  by  mea.su  rinj 
the  frequency  shift  and  change  in  Q  caused  by  introducing  the  plasma  into 
the  microwave  cavity.  The  equation  used  to  find  the  electron  density  and 
collision  frequency  are  based  on  a  perturbation  solution  that  assumes  that 
the  fields  in  the  presence  of  the  plasma  are  not  very  different  from  those 
of  the  unperturbed  cavity.  When  the  plasma  density  is  low  enough  for  tin: 
approximation  to  be  valid,  the  discharge  current,  which  is  proportional  to 
electron  density,  and  the  resonant  frequency  shift  have  an  approximately 
linear  relationship.  Figure  23  shows  the  detuning  of  the  microwave  cavity 
as  a  function  of  the  discharge  current  for  fixed  hydrogen  pressure  and 
fixed  magnetic  field.  The  relationship  is  linear,  so  the  perturbation  tech¬ 
nique  will  be  viseful  in  establishing  the  operating  point  of  the  plasma. 
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Additional  diagnostic  data  are  given  in  associated  Report  3  and  show  the 
effect  of  varying  magnetic  field.  Current  voltage  characteristics  for 
the  discharge  at  various  values  of  magnetic  field  are  also  given  there. 

The  discharge  was  found  to  be  quite  stable  and  the  plasma  was 
easily  controlled  by  varying  the  discharge  current.  It  was  also  pos¬ 
sible  to  reset  the  current  so  that  the  cavity  detuning  was  within  one  or 
two  megacycles  of  the  previously  obtained  values. 

The  electron  beam  operated  at  5  kv  with  a  microperVeance  of 
approximately  unity  and  operated  normally  in  a  solenoid  magnetic  field 

_3 

with  or  without  hydrogen  pressure  of  up  to  10  mm  Hg  . 

The  gain  of  the  amplifier  system  operating  as  a  klystron  has  been 
calculated  to  serve  as  a  starting  point  for  the  investigation.  The  beam 
diameter  (0.180  in.)  was  very  much  smaller  than  the  diameter  of  the 
drift  tube  (0.650  in.).  Thus,  the  coupling  to  the  beam  was  expected  to 
be  low  . 

5 

The  gap  coupling  coefficient  was  calculated  using  the  expression, 


!o  ^b)  "  A  1 

ij  (va)  d2 


F(z)  cos  [}^z  dz  , 


F(z)  -  Ez/Em-  is  the  normalized  gap  field.  All  other  symbols  have  their 
usual  meaning.  The  field  distribution  was  measured  experimentally  and 
the  integral  evaluated  graphically,  where 

d  =  Drift  space  length  =  5.08  cm  , 

b  =  Beam  radius  =  0.090  in.  =  2.28  x  10  ^  m  , 
a  =  Radius  of  drift  tube  =  0.325  in.  =  8.25  x  10  ^  m 
Hence,  =  1.02  x  10  ^  . 


-49- 


pression. 


The  gain  of  a  two  cavity  klystron  amplifier  is  given  by  the  ex- 
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F  =  Reduction  factor  =  1 
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R  =  — —  the  d-c  beam  impedance  =  20  KS2  , 
° 


G—  =  +  — —  =  0.44  x  10  5  mho  , 

T  Rsh  2Ro 


R  ^  “  Shunt  resistance  of  the  cavity  =  125  K  , 


■  3.06  x  lo”  c/ s  , 


(G)  db  -  10  log  G  =  -24  db 
The  klystron  gain  measured  experimentally  was  -30  db 
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VII.  CONCLUSIONS 


It  can  be  concluded  that  a  test  vehicle,  having  a  plasma  generating 
chamber  separate  from  an  electron  gun  chamber,  can  be  constructed  to 
yield  more  than  100:1  pressure  ratios  to  allow  dense  plasma,  on  the 
one  hand,  and  a  low  pressure  environment  for  the  electron  gun,  on  the 
other  hand,  and  that  high-power  convergent  gun  electron  beams  operate 
weil  under  these  conditions. 

It  can  also  be  concluded  that  a  flexible  subassembly  can  be  con¬ 
structed  in  the  plasma  chamber  to  allow  the  generation  of  a  plasma  and 
input  and  output  direct  microwave  coupling  to  a  plasma. 

It  can  further  be  concluded  that  a  beam  scanner  can  be  constructed 
to  study  an  electron  beam,  sent  through  a  gas  atmosphere,  for  fluctua¬ 
tions,  as  well  as  to  study  the  microwave  properties  of  such  an  electron 
beam  after  efficient  interaction  has  taken  place, 

Finally,  it  can  be  concluded  that  the  test  vehicle  developed  can 
be  used  effectively  to  study  the  interaction  of  a  pulsed  electron  beam  of 
a  few  kilowatts  with  a  controlled  plasma  of  a  double  hot-cathode  Penning 
configuration  at  3,000  Mc/s  operating  frequency. 
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ANALYSIS.  OF  THE  INTERACTION 
OF  AN  ELECTRON  BEAM  WITH  A  PLASMA 

H.  Marantz 
SUMMARY 

The  theoretical  anaLysis  undertaken  in  this  report  has  two  ob¬ 
jectives.  First,  to  understand  the  physical  processes  of  plasma  genera¬ 
tion  in  a  low-pressure  Penning  discharge.  Second,  to  obtain  quantitative 
information  on  the  factors  which,  independent  of  external  coupling,  limit 
the  gain  of  the  beam-plasma  interaction. 

In  the  low-pressure  double -Penning  discharge,  particle  losses 
from  the  system  are  due  more  to  the  effects  of  free  -  streaming  and  elec¬ 
tric-field  acceleration  than  to  diffusion.  This  is  the  case  whenever  the 
axial  moan-free-path  is  comparable  to  the  length  of  the  system.  The  dis¬ 
charge  characteristics  predicted  by  this  model  are  much  closer  to  the 
experimental  data  than  those  predicted  by  the  high-pressure  model  with 
axial  diffusion. 

The  beam -plasma  interaction  is  handled  both  qualitatively  and 
quantitatively.  The  explanation  of  the  interaction  mechanism  is  developed, 
first  from  the  principles  of  particle  trapping  and  energy  conservation  and, 
then,  in  terms  of  the  familiar  microwave  device,  the  multi -cavity  klystron. 
The  interaction  is  treated  quantitatively  by  using  appropriate  integrals 
of  the  Boltzmann  equation  to  describe  the  microwave  characteristics.  Thi 
information,  together  with  Maxwell's  equations,  form  a  closed  set  of 
equations  from  which  the  fields  are  obtained.  The  analytical  approach  for 
solving  the  eigenvalue  problem  is  formidable  for  all  but  the  simplest 


approximations .  We  have  closed-form  solutions  for  these, and  numerical 
results  were  obtained  for  the  remaining  cases  with  a  digital  computer. 

Analysis  of  the  interaction,  using  a  simple  one-dimensional  model, 
showed  that  extremely  high  (>100  db/cm)  gain  was  possible.  The  results 
include  the  effects  of  finite-transverse-inhomogeneous  beam-plasma 
cross  section,  magnetic  field,  plasma  collision,  and  plasma-electron 
temperature.  These  conditions  are  all  present  in  a  real  plasma  and  serve 
to  limit  the  gain  to  about  10  db/cm  at  S-band  frequencies.  Although  this 
is  an  order  of  magnitude  less  than  that  predicted  by  the  simplest  model, 
it  is  still  higher  than  existing  devices  at  these  frequencies. 
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PLASMA  DIAGNOSTICS 


D.  M.  Kerr 
SUMMARY 

Two  singly  re-entrant  cavities  of  the  type  ordinarily  used  in 
klystron  amplifiers  were  designed  and  tested  for  use  in  the  electron- 
beam  hydrogen-plasma  system.  They  were  designed  primarily  for 
use  as  the  input  and  output  coupling  structures  of  the  beam-plasma 
amplifier,  but  they  also  are  used  to  investigate  the  properties  of  the 
plasma  discharge,  in  particular,  to  determine  the  electron  density  and 
collision  frequency. 

The  design  and  cold  tests  of  the  cavities  are  described,  and 
comparisons  are  made  between  theoretical  and  experimental  values. 

A  theoretical  discussion  of  cavity  Q  in  the  presence  of  coupling  losses 
is  included,  since  the  measurement  of  cavity  Q  is  important  in  diag¬ 
nosing  the  plasma.  A  perturbation  technique  was  used  to  determine  the 
quantity,  p^Rg^/Qo  ,  which  is  important  in  considering  coupling  to  the 
system. 

The  general  theory  of  the  cavity  perturbation  is  given  and  equa¬ 
tions  are  obtained  relating  the  shifts  in  resonant  frequency  and  Q  to  the 
average  electron  density  and  collision  frequency.  The  possible  dis¬ 
advantages  of  this  type  of  measurement  are  discussed  briefly  and  then 
preliminary  data  is  presented.  A  linear  relationship  between  resonant 
frequency  shift  and  discharge  current  is  found,  as  predicted,  and,  in 
the  range  of  discharge  currents  expected  in  the  beam-plasma  amplifier, 
this  type  of  diagnostic  technique  appears  to  be  successful. 
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ION-INDUCED  OSCILLATIONS  OF  ELECTRON  BEAMS 


A.  S.  Gilmour,  Jr, 

SUMMARY 

In  this  study,  the  ion-induced  oscillations  of  a  pulsed,  conver¬ 
gent  flow,  Brillouin  beam  were  measured  in  detail  for  pressures  of 
“8  -  3 

hydrogen  gas  from  10  Torr  to  10  Torr  .  The  beam  passing  through 
the  hydrogen  was  .2.0  in.  in  diameter  and  had  ,40  a  current  for  5  kv  . 

The  beam  dynamics,  as  they  depended  on  time  during  the  pulse,  were 
measured,  using  a  beam  scanner  that  was  able  to  measure  the  beam 
current-density  at  any  transverse  location  and  over  a  considerable  length 
of  the  beam.  Associated  test  equipment  allowed  sampling  of  any  desired 
portion  of  the  time  pulse.  The  experimental  results,  showing  definite 
transverse  and  radial  ion-induced  electron-beam  oscillations,  are  given 
in  detail  and  explained.  The  transverse  oscillations  occurred  at  the  lower 
pressures  and  increased  with  time  during  the  pulse.  The  radial  oscilla¬ 
tions  occurred  at  the  larger  pressures  and  also  increased  with  time  during 
the  pulse.  It  was  concluded  that  the  feedback  mechanisms  for  these  two 
forms  of  oscillation  are  the  charged  particles  moving  from  the  drift 
region  back  toward  the  gun  region  and  affecting  the  beam-formation  con¬ 
ditions  there. 

Earlier  work  by  Mihran  and  Agdur  is  discussed  and  extended,  and 
the  theoretical  work  of  Peterson  and  Puthoff  is  extended  to  explain  the 
oscillation  build-up. 

Ways  of  controlling  or  eliminating  these  oscillations  are  deduced 
and  explained. 
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ABSTRACT 


A  detailed  analysis  is  made  of  the  gain  characteristic,  of  the 
beam-plasma  amplifier.  The  first  section  is  a  study  of  the  general 
specification  of  a  plasma  and  of  the  d-c  plasma  discharge  in  partic¬ 
ular.  It  is  found  that  the  analysis  of  the  low-pressure  Penning  dis¬ 
charge  must  take  particle  streaming  under  consideration  to  account 
for  the  discharge  current  and  the  percentage  of  ionization. 

The  second  section  analyzes  several  mechanisms  for  broad¬ 
ening  the  gain  characteristic  of  an  infinite  beam  plasma,  it  is  con¬ 
cluded  that  the  beam-plasma  inhomogeneities  account  for  the  great¬ 
est  change  in  the  gam  characteristic.  Corrections  for  plasma  col¬ 
lisions  and  temperature  can  be  neglected  for  the  experimental  values 
expected. 
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or 


whether  the  parameter  is  a  function  of  time.  Double  subscripts  indicate 
both  of  these  factors.  The  first  subscript  is  used  for  time  dependence, 
the  second  for  particle  type. 


I.  INTRODUCTION 


Interest  in  plasmas  started  in  1902  with  Heaviside's  postulation 
of  an  ionized  layer  above  the  earth.  Since  then,  two  major  fields  of 
interest  have  developed:  first,  the  use  of  a  plasma  as  an  energy  con¬ 
verter,  e.g.  ,  thermonuclear  machines  and  ion-propulsion  guns;  and 
second,  the  use  of  a  plasma  as  a  wave -propagating  medium,  e.  g.  , 
magneto -ionic  theory  and  the  production  of  shock  and  electromagnetic 
waves  in  a  plasma. 

We  are  concerned  with  using  a  plasma  to  amplify  and  generate 
communication  signals  at  about  3  Gc/s  .  The  possible  advantages  of 
using  a  plasma  are  fourfold:  (1)  it  allows  propagation  of  power  at  high 
frequencies,  where  metallic  slew-wave  structures  are  difficult  to  con¬ 
struct;  (Z)  it  has  strong  interaction  with  an  electron  beam,  yielding 
high  electronic  gain  per  centimeter  of  drift  length;  (3)  electronic  tun¬ 
ing  of  the  plasma  is  possible;  and  (4)  the  nonlinear  properties  of  the 
plasma  could  be  used  in  frequency  multiplying  applications. 

Pierce  was  one  of  the  first  to  think  of  using  a  plasma  in  a  prac¬ 
tical  device  for  microwave  amplification.  The  basic  interaction  theory 

Z  3.4 

was  developed  in  papers  by  Bohm  and  Gross,  Haeff,  and  Pierce. 

5  6 

More  recently,  the  analysis  has  been  extended  by  Trivelpiece,  Paik, 

7  8 

Chorney,  and  Bers  to  include  the  effects  of  finite  geometry.  The 

theory  of  the  double  stream  instability  for  a  nonzero  temperature  plasma 

9  10 

has  been  discussed  by  E.  A.  Jackson,  J.  D.  Jackson,  and  Kellogg 
11 

and  JLimohn. 
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In  this  study,  a.  Penning  discharge  device  and  a  microwavetube 
electron  gun  were  placed  back  to  back  (see  Figures  1  and  2).  This 
discharge,  by  using  a  confining  axial  magnetic  field,  generated  high 
plasma  densities  at  relatively  low  gas  pressures.  The  electron  gun 
used  semiconfined  focusing  to  send  a  5-kv,  .  4-a  (10-p.s  pulsed)  beam 
through  the  plasma  region.  Microwave  signals  were  directly  coupled 
to  and  from  the  beam-plasma  region  by  using  singly  re-entrant  klystron 


cavities. 


II.  CHARACTERISTICS  OF  A  PLASMA 


A.  GENERAL  CHARACTERISTICS 


A  plasma  is  an  ionized  gas  in  which  the  positive  particle  den¬ 
sity  is  approximately  equal  to  the  negative  particle  density.  Plasmas 
occur  in  nature  (sun's  corona,  earth's  ionosphere)  and  in  man-made 
devices  (arc  discharges  and  flame  of  Bunsen  burner).  There  are  four 
parameters  that  determine  the  assumptions  to  be  made  in  our  analysis. 
First,  the  Debye  length  h  must  be  considered.  This  is  the  shielding 
distance  for  the  potential  resulting  from  a  single  charged  particle. 
Therefore,  even  though  we  assume  equal  numbers  of  positive  and  nega¬ 
tive  charges,  potentials  still  exist  for  distances  of  the  order  of  the 
Debye  length,  which  is  given  by 


h2  = 

l 


kT. 

1 


(4tt  m  ql*) 


(2.  1) 


Second,  the  average  distance  between  charged  particles  d  is 
d  =  (2n)-1/3 


(2.2) 


Third,  the.  average  critical-impact  parameters  for  the  collision  of  two 
charged  particles  P^  is  defined  as  that  impact  parameter  which  re¬ 
sults  in  a  scattering  angle  of  tt/2  .  For  the  coulomb  potential,  we  get 


P  . 

oi 


3kT- 


(2.  3) 
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The  fourth  parameter  is  the  average  mean-free  path  between  strong 


charged-particle  collisions  •  If  we  use  a  coulomb  potential,  we  get 


2  -1 

=  (im.P  .) 

1  oi 


When  the  Debye  length  is  much  smaller  than  a  typical  dimen¬ 
sion  of  our  geometry,  we  know  that  in  the  absence  of  strong  external 
fields,  the  plasma  volume  can  be  assumed  to  be  electrically  neutral. 

When  Pq  «  d,  the  close  collisions  are  strictly  two-particle  colli¬ 
sions.  This  justifies  using  the  coulomb  or  modified  coulomb  potential 
in  calculating  the  effects  of  plasma  collisions.  If  Pq  >  d,  the  concept 
of  a  distinct  two -particle  collision  must  be  replaced  by  that  of  collective 
interaction,  and  the  effects  of  collisions  cannot  be  easily  isolated.  For 
d  «  X  ,  the  effective  mean-free  path  (distance  between  deflections 
of  ninety  degrees  or  more)  is  long.  The  expansion  of  the  collision 

term  in  terms  of  small-velocity  changes  leads  to  the  Fokker  -  Planck 
12 

Equation.  In  order  to  determine  the  d-c  properties  of  the  discharge, 
we  must  know  the  particle  density  and  the  root-mean-square  particle 
velocity.  For  the  microwave  properties,  we  must  be  able  to  calculate 
the  r-f  current  and  space-charge  density.  The  approach  used  to  get 
these  quantities  depends  on  the  time  scale  for  which  the  information 
must  be  valid. 

i)  During  a  collision,  y/e  use  the  BBKGY  relations.  These 
equations  are  derived  by  taking  successive  moments  of  Liouville's  equa¬ 
tion  with  respect  to  the  velocity  of  each  particle.  For  a  plasma  with 
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n  negative  particles,  this  involves  2n  coupled  differential  equa¬ 
tions  . 

ii)  Between  collisions,  we  use  Boltzmann's  equation;  this  in¬ 
volves  only  two  coupled  differential  equations. 

iii)  After  several  collisions,  we  use  the  moment  equations  that 
are  derived  from  Boltzmann's  equations.  There  are  four  of  them, 
but  we  have  three  fewer  independent  variables  than  in  item  ii.  Typi¬ 
cally,  the  times  for  this  experiment  are: 

- 15 

10  s  for  collision  time, 

10  1  s  for  the  txm  e  between  collisions, 

_Q 

10  s  for  the  r-f  time  of  interest, 

-  5 

10  s  for  the  d-c  time  of  interest  . 

We  shall  therefore  use  the  Boltzmann  transport  equations  to  describe 
the  d-c  discharge  and  the  Boltzmann  equation  to  describe  the  micro- 
wave  phenomena,  These  equations  are  derived  in  Appendix  A. 

B.  THEORY  OF  THE  PLASMA  DISCHARGE 

A  low-voltage,  hydrogen-plasma  (thermionic)  discharge  was 
used,  the  geometry  of  which  is  shown  in  Figure  2.  The  cathode  elec¬ 
trons  were  accelerated  by  the  cathode-anode  potential,  so  that  they  had 
an  energy  greater  than  the  ionization  potential  when  they  passed  through 
the  center.  Since  these  electrons  collide  with  the  neutral  molecules, 
they  never  return  to  the  cathodes  and  must  eventually  reach  the  anode. 
Some  -f  these  give  up  their  energy  in  ionizing  collisions  before  being 
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collected  by  the  anode.  The  plasma  of  electrons  and  Hj  thus  created 
will  "diffuse"  away  from  the  point  of  ionization,  the  ions  going  to 
the  cathode,  the  electrons  to  the  anode.  The  term  "diffuse"  is 
used  because  it  is  not  diffusion  in  the  usual  sense  that  governs  the 
longitudinal  flow  of  particles  in  our  system.  When  the  system  pres¬ 
sure  is  low  enough  that  the  particle  mean-free  path  is  of  the  same 
order  of  magnitude  as  a  typical  dimension  of  the  confining  structure, 

particle  streaming  and  electric -field  acceleration  must  be  taken  into 
13 

account.  Simon  has  used  these  terms  to  explain  anomalous  diffusion 
in  arc  plasmas.  The  equations  governing  the  discharge  are  (see  Ap¬ 
pendix  A): 


V+  =  V  v  •  -£+•  v 


_  %nn+  eEn 

’  n  +  u  •  E  n  -  - — -  •  - — ~ 

+  5  +  .  -  +]  l  m+v  + 


d+n  =1-0 


•  j^-D  _  •  Vn_+ji  ■  En_  *•- 


V-Hn-  ,  CE! 


y  •  e^E  =  e(n+  -  n_) 


(2.7) 


In  equations  (2.  5)  and  (2.  6),  -v  ^n_yL  represents  the  free 
streaming  effects  and  e  E  |  m  v  |  j  the  electric -field  acceleration 
effects.  Note  that  these  effects  are  not  important  in  the  transverse 
direction  because  of  the  strong  magnetic  field  required  to  confine  the 
electron  beam.  Therefore,  we  do  not  include  the  transverse  compo  - 
nents  of  these  terms  in  Equations  (2,  5)  and  (2.6).  Conversely,  because 
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of  the  long  axial  mean  -free  path  at  low  pressures,  the  axial  diffusion  terms 
do  not  apply.  If  we  calculate  the  highest  possible  electric  field  that 
could  be  present  in  the  plasma  in  the  absence  of  oscillations,  we  find 

E  W  iO4  v/m 
max  ' 

Using  Gauss's  law,  we  calculate  the  maximum  deviation  from  a  neu¬ 
tral  plasma,  which  is 
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Multiply  Equation  (2.9)  by  ai+v+ j  j  ,  Equation  (2.  10)  by  m  v  j  j  ,  and 
add 


V  = 


Vlll1!  +  m-V- 
m+V+  |  |  +m-' 


m , v .  ,  ,D  ,  , 

.  +  m  v  , 

i  D  .I 

+  +1 1  +1 

1  -  -  1 

i  - 1 1 

"VMl 

+  m-V-|l 

r  2 

m  .  v  . 


+  +  i  i 


+  m  v 


m+v+ll  +  m-v.|| 


By  setting  the  left-hand  side  of  Equation  (2.  9)  equal  to  zero, 

calculating  I  1^  ,  and  assuming  appropriate  boundary  conditions ,  we 

could  solve  for  the  exact  steady-state  solutions.  We  can  obtain  much 

information,  however,  using  much  simpler  calculations  and  heuristic 

reasoning.  Most  of  the  ionizations  take  place  in  the  volume  connecting 

the  two  plasma  cathodes  The  originally  empty  inner  cylinder  becomes 

filled  with  plasma  by  transverse  diffusion.  If  we  neglected  the  axial 

particle  motion,  we  could  assume  that  the  particle  density  was  constant 

across  the  inner  cylinder  diameter.  In  fact,  Chorney  etal.  ,  using 

a  pair  of  annular  cathodes  under  much  the  same  conditions  as  ours, 

have  observed  that  the  plasma  light  intensity  was  lowest  in  the  center. 

This  does  not  necessarily  mean  that  the  inner  cylinder  has  a  lower 

electron  density,  but  that  the  ratio  of  fast  to  slow  electrons  is  smaller 
14 

there.  Chorney  later  noted  that  a  discharge  instability  occurred  at 

high  magnetic  fields.  This  effect  has  been  analyzed  by  Kadomtser 
15 

etal.  They  found  that  the  instability  was  associated  with  the  growth 

of  a  low-frequency  6-dependent  mode.  In  Chorney' s  experiment,  the 
plasma  was  enclosed  by  glass  walls;  therefore  the  excitation  of  an 
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axially  unsymmetric  mode  was  much  more  probable  than  in  our  own 
case.  Therefore  we  do  not  expect  a  critical  magnetic  field  that  would 
be  low  enough  to  interfere  with  good  electron  beam  focusing.  Figure  3 
shows  the  particle  density  cross  section  calculated  using  a  one-dimen¬ 
sional  model.  Also  shown  there  are  some  analytical  expressions  that 
could  be  used  to  approximate  this  distribution  profile.  For  d-c  cal¬ 
culations  the  approximation  using  the  series 


—  =  V  C.J  (b.i 
m  ZL  i  o'  i 


where  J  is  the  zero  order  Bessel  function  of  the  first  kind,  is  very 
useful.  Setting  I  =  I  =  0  in  Equation  (2.  11),  we  can  calculate  the 
amount  of  discharge  current.  Solving  Equation  (2.  11)  by  separation  of 
variables  and  neglecting  all  but  axially  symmetric  modes,  we  find  that 
the  terms  in  the  solution  are  of  the  form, 


n  =  n  .  e  J  (b.r) 
oi  o'  1 


(2.  13) 


where 

X 

The  n^j  are  determined  by  the  initial  (steady-state)  configuration.  The 
current  is  given  by 


m+v+ 1|  +  m-v- 1  | 


Kvdn 

L[m+v+| 


2  2  I 

HI+m-V-||. 


•  (2.  14) 


I  =  2  |  n  dv  1  -  —  1 

e\(b2) 

(2.  15) 

■>  l  eJ 
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FICiURK  3.  Plasma  Cross  St*rliom>. 


b2  =  — - -  .  (2.16) 

Z”o. 
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Implied  in  this  derivation  is  the  assumption  that  the  presence  of  ion¬ 
ization  processes  will  not  significantly  change  the  rate  at  which  parti¬ 
cles  leave  the  system.  There  will  not  be  large  gradients  in  the  spatial 
distribution  of  neutral  particles,  so  that  the  minimum  pressure  that 
will  supply  all  the  charged  particles  is  given  by 

p  .  =  n  kT..  x  .75  x  10  2  mm  Hg  (cgs)  .  (2.17) 

The  actual  pressure  is  calculated,  knowing  Pq  ,  the  probability  of  ioni¬ 
zation  per  centimeter  per  millimeter  of  mercury,  the  mean  electron 
velocity,  and  the  time  it  takes  a  fast  electron  to  leave  the  system: 


For  work  on  hot  plasmas,  it  is  desirable  to  know  the  ratio  between 
fast  and  slow  electrons,  so  that  the  electron  velocity  distribution  can  be 
approximated.  First,  we  calculate  a  transverse  electron  drift  velocity 
by  considering  the  initial  position  of  the  last  set  of  electrons  to  reach 
the  anode  after  1/X  seconds.  The  ratio  of  the  time  the  electron  spends 
as  a  fast  electron  to  that  as  a  slow  electron  gives  us  the  ratio  of  num¬ 
ber  of  fast  to  number  of  slow  electrons  in  the  volume. 
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Using  values  expected,  in  our  experiment,  we  should  be  able  to 
predict,  to  within  an  order  cf  magnitude,  discharge  current,  percent¬ 
age  of  ionization,  and  the  ratio  of  faBt  to  slow  electrons. 

Using  the  discharge  geometry  of  Figure  2  and  the  cross-section 
of  Figure  3,  we  compute  a  typical  set  of  discharge  parameters.  Given 


VD 

15  v 

u  = 

PF 

2 it  x  3  kmr/s 

0)  = 
ce 

2 ir  x  2  25kmr/ 

= 

7 

10  kmr/s 

r  - 

o 

.75  cm 

I* 

10  cm 

tn 

10  3  K 

b2 

2  x  104  m'2 

we  calculate 
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III.  THEORY  OF  THE  INTERACTION 


A.  PHYSICAL  EXPLANATION 


A  plasma  can  support  many  different  types  of  waves.  These  are 
usually  described  in  two  ways:  (1)  whether  the  electrons  or  the  ions 
or  both  provide  the  particle  reaction  that  comprises  the  wave,  (2)  the 
type  of  electric  or  magnetic  field  that  can  be  associated  with  the  wave. 

We  are  concerned  with  the  electrostatic  waves,  which  have  their  elec¬ 
tric  fields  parallel  tc  the  direction  of  propagation,  and  the  electro¬ 
magnetic  waves  with  their  electric  fields  perpendicular  to  the  direction 
of  propagation.  Both  are  also  high-frequency  waves;  i.  e.  ,  the  elec¬ 
trons  oscillate  back  and  forth  in  a  background  of  stationary  positive  ions. 

The  amplification  mechanism  is  explained  by  following  an  elec¬ 
trostatic  wave  as  it  propagates  through  the  plasma.  The  wave  travels 
at  phase  velocity  vq  with  a  finite  amplitude;  therefore,  the  wave 
trough  can  trap  those  electrons  whose  streaming  velocities  are  in  a 
range  v  +  A  <  v  <  v  +  A.  Once  trapped,  the  electrons  move  with 
the  wave  trough  at  velocity  vq  .  Now  consider  the  electron  distribution 

function  f.  at  velocities  v  +  A  .  Let  6  =  f.(v  +  A)  -f.(v  -  A). 

1  O  —  1  o  1  o 

If  6  ?■  0  ,  there  is  a  net  transfer  of  energy  between  the  wave  and  the 
particles,  since  the  total  energy  of  the  system  is  conserved.  For  small 
signals,  the  net  energy  transferred  is  given  by 


AE 


m.  f .  (v  )  v  A 
i  l  o  o 


8f(v) 

3v 


v  = 


V 
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Therefore,  if  the  slope  of  L(v)  is  positive  at  vq  ,  the  transfer  of 
energy  is  to  the  wave,  and  the  trough  grows  in  amplitude,  traps  more 
particles,  etc.  until  some  nonlinearity  levels  off  the  growth.  Typi¬ 
cally,  this  could  be  due  to  particle  overtaking;  i.  e.  ,  the  r-f  electron 
velocity  becomes  greater  than  the  d-c  electron  velocity.  If  the  slope 
is  negative  at  v  ,  the  energy  is  transferred  to  the  particles,  and  the 
wave  amplitude  decays  until  it  reaches  the  noise  level  associated  with 
the  system.  This  is  a  form  of  Landau  damping. 

In  microwave  terms,  the  plasma  is  analogous  to  a  sequence  of 
uncoupled  cavities,  all  resonant  at  the  same  frequency.  When  a  beam 
is  passed  through  these  cavities,  any  beam  disturbance  that  has  a 
Fourier  frequency  component  near  the  resonant  frequency  of  the  cavity 
will  excite  an  oscillation  in  the  cavity,  which  will  simultaneously  modu¬ 
late  the  beam.  In  a  multicavity  klystron  for  example,  the  first  cavity 
modulates  the  velocities  of  the  beam  electrons  with  an  external  micro- 
wave  signal.  The  beam  then  travels  through  a  drift  space  where  the 
electrons  become  density  modulated.  When  the  beam  arrives  at  the 
second  cavity,  it  excites  an  oscillation  that  simultaneously  remodulates 
the  beam.  Each  subsequent  cavity  save  the  last  one  impresses  a 
stronger  modulation  on  the  beam,  thus  converting  an  increasing  amount 
of  the  kinetic  energy  of  the  beam  to  r-f  energy.  In  the  beam-plasma 
amplifier,  the  cavities  are  distributed  rather  than  discrete;  therefore 
the  signal  gain  is  exponential  as  in  a  traveling -wave  tube.  Unlike  the 
usual  TWT  interaction  structure,  however,  the  plasma  does  not  nec¬ 
essarily  support  a  propagating  wave  by  itself. 
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These  are  two  points  to  be  made  in  interpreting  the  complex 
propagation  constant.  First,  we  can  expect  gain  only  when  the  phase 
velocity  of  the  waves  is  less  than  that  of  the  d-c  velocity.  Second, 
we  are  primarily  interested  in  the  forward-wave  interaction;  i.e,  ,  the 
velocity  of  energy  transport  is  in  the  same  direction  as  the  phase  veloc¬ 
ity  of  the  wave.  The  backward-wave  interaction  is  also  present  for 
certain  frequencies,  but  for  high-frequency  amplifiers,  the  forward- 
wave  interaction  has  greater  bandwidth  and  stability.  The  gain  char¬ 
acteristics  (Figures  5-7,  9-13)  usually  show  both  forward  and  backward 

2  ~ 

waves,  especially  for  1  <  a  <  (1+d^)  ;  therefore,  v^.  must  pass 

through  zero.  As  a  result  we  get  a  zero  in  the  actual  gain  at  the 

boundary  between  these  two  types  of  waves.  We  shall  set  v-g  =  , 

although  this  is  not  strictly  valid  in  a  region  where  the  propagation 

17 

constant  is  complex  and  varies  rapidly. 

These  explanations  are  most  applicable  and  easiest  to  use  for 
the  ideal  case  of  infinite,  homogeneous  beam-plasma  cross  sections. 
They  do  not  include  the  effects  of  finite  geometry,  magnetic  field,  in¬ 
homogeneous  cross  sections,  collisions,  and  nonzero  temperature . 
These  effects,  which  are  present  in  a  real  plasma,  are  treated  in  the 
following  section. 

B.  MATHEMATICAL  DESCRIPTION 

i.  General  Equations 

We  use  three  sets  of  equations  to  describe  the  interaction.  The 
first  set  describes  the  wave  phenomena  as  postulated  by  Maxwell; 
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xE  =  -pd  H  , 

—  *0  +  — 

(3.  1) 

xH  =  e  d  E  +  J  i 

(3.2) 

—  o+  — 

•  Co  E  =  P 

(3.3) 

•  p  H  =  0 

(3.4) 

The  second  set,  postulated  by  Boltzmann,  describes  the  properties  of 
the  media  (see  Appendix  A): 

q. 

d,f.  +  (v  .  •  D)f.  +  -i-  f(E  +  V.  x  B)  -  v  V  If.  =  0  .  (3.  5) 

+  i  '—i  l  m.  L—  ~1  1J  1 

The  subscript  i  refers  to  the  group  of  particles  of  type  i  . 

The  third  set  of  equations  couples  the  two  previous  ones: 

-  =  Zqi/VifidV>  ’  (3’6) 

i 

p  =  Zqi/fi  dvi  •  (3'7) 

i 

These  equations  will  be  solved  subject  to  the  following  boundary  condi¬ 
tions.  The  boundary  is  a  perfectly  conducting  cylinder,  infinite  in 
length,  and  solutions  for  the  axial  electric  field  must  be  square  inte- 
grable  over  any  unit  cross-sectional  area.  This  means  that 
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We  are  interested  in  finding  the  steady-state  axially  symmetric  modes. 
These  modes  have  a  nonzero  axial,  electric  field  at  the  center  of  the 
cylinder,  thus  ensuring  a  strong  interaction  between  the  beam  and  the 
plasma.  Therefore,  we  look  for  solutions  that  have  a  space-co-ordinate 
dependence  of  the  form  R(r)  exp  j(ut  -  j3z)  .  This  simplifies  the  fol¬ 
lowing  partial  derivatives  to 

d+  =  ^  dz  =  “jP  de  =  0  • 


The  media  equations  introduce  the  nonlinear  term, 


x  B) 


Though  the  solution  with  nonlinearity  included  is  a  problem  of  great 
interest,  we  consider  the  linearized  equations  to  keep  the  mathematics 
simple.  To  linearize  the  equations,  we  expand  all  the  variables,  except 
time,  in  the  following  manner: 


A(x,  v.  t)  =  A  (x  ,  v  )  +  \A(x  ,  v  ,  t) 

—  o  — o  — o  o  o 


(3.6) 


where  \  is  an  indication  of  the  order  of  nonlinearity.  For  the  lin¬ 
earized  solutions  we  require  that 


Put  Equation  (3.3)  into  the  three  sets  of  equations,  then  collect  terms 
by  their  power  in  \  .  By  making  \  arbitrarily  small,  we  can  justify 
setting  each  coefficient  of-  k1  equal  to  zero,  thereby  getting  a  zero- 
order  and  a  first-order  set  of  equations,  which  are 


-Z0- 


(3.  23) 


Equations  (3  10)  through  (3.  to)  form  a  nonlinear  time  -  independent,  set. 
Eqea’irns  '3.  I7)  *hrcugh  (3.23*  form  a  linear  time -dependent  set.  As¬ 
sume  that  ^  0  This  ij  a  reasonable  assumption  as  long  as  the 

17  -  3 

charged  parvcle  density  ^  0  {This  requires  that  p  >  10  m  ). 
As;i,me  that  we  have  solved  the  time -independent  set  of  equations.  We 
usual’v  these  t he  aero  erde-  distribution  functions  to  be  Maxwellian. 
This  is  reasonable  because  the  electron-beam  duty  cycle,  i.o  ,  ratio 
of  on  ic  off  time.  is  very  low.  and  Boltzmann's  H  theorem  states 
that  *he  equilibrium  distribution  function  is  Maxwellian.  To  solve  the 
time  -  dependent  sc-t.  we  express  f  in  terms  of  the  variables,  , 

H  .  v.  .  and  then  evaluate  either  J  l  or  pt  by  performing  the  appro¬ 
priate  ir.teizrais  ever  velocity  space  We  put  these  results  either  into 
the  wave  equation  or  into  Poisson's  equation  and  solve  for  the  sclf-cor- 
sistcnt  t  r.ontr  :v:a!)  fields  This  gives  us  a  differential  equation  whose 
solution;  are  eigenfunctions,  from  which  we  get  the  dependence  of  the 
fieid  on  the  spatial  co-ordinates  We  also  get  the  longitudinal  propaga¬ 
tion  constant  p  as  a  function  of  frequency. 


2  Effects  of  Geometry  and  Magnetic  Field 


We  take  f  to  be  a  normalized  shifted  Maxwellian  distribution 
oi 

m  the  zero  order  longitudinal  velocity  v.  .  This  gives 

.2 


.  .  .  -i/i  '-i/i  -  '•<,,> 

'of'.1  ‘  -  VT  c - X- 


(3. 24) 
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However,  we  will  anticipate  the  results  of  Section  3b  which  follows, 

that,  a  nonzero  plasma  temperature  dees  not  strongly  affect  the  gain 

characteristic  ’-hen  the  geometry  is  finite,  and  consider  the  case  where 

v  approaches  zcrc:  then  f  .  becomes  the  Dirac  delta  function  de- 
Tj  01 

tmed  bv 


oc 

1  =  /  foi(Vj)  dVi 

“CO 


(3.  25) 


This  is  easily  seen  by  making  the  change  of  variables, 


v  .  -  v  . 

U  =  — - —  (3.26) 

v  T. 

i 

ard  taking  the  integral  of  the  limit.  Since  f  is  defined  by  an  inte- 

oi 

gral  over  velocity  space,  we  must  use  the  velocity  moments  of  the 

Boltzmann  equation  instead  of  the  equation  itself  to  describe  the  beam- 

18 

plasma  media.  Following  the  derivation  given  by  Spitzer,  we  get  the 
following  equations: 


dj.  v  + 

+  —  li 

(v.-v)v  +(v  *v)v.=F. 

' - Ol  '  —  ll  —  ll  —  Ol  1 

(3.  27) 

F.  =  N.  1 
1  1  1 

E  ,  +  v  ,  x  B  +  v  .  x  B  .  1 
—  1  — ll  — o  — Ol  —  1J 

(3.  28) 

d  N  . .  + 

+  li 

V  •  (JN.v  +N  v,)  =  0 

'  01  — ll  ll  Ol 

(3.29) 

(<H 

11 

(N  .  q.  v  +  N  q.  v  ,) 

oi  ^i  —  li  li  n  —  oi 

(3.  30) 

1 

q.  N  . . 
n  li 

(3.  31) 

i 


-  23- 


Equations  (3.  27)  -  (3.  31),  the  wave  equations,  and  the  boundary  condi¬ 
tions  give  enough  information  to  solve  both  for  the  propagation  constant 
as  a  function  of  frequency,  and  also  for  the  r-f  fields.  In  the  general 
case,  however,  this  involves  solution  of  two  coupled,  linear,  fourth- 
order  differential  equations.  There  are  two  types  of  approximations 
that  reduce  the  order  of  the  differential  equation  from  four  to  two. 

The  mode  decoupling  approximation  is  used  when  the  waves  of  inter¬ 
est  would  propagate  in  the  absence  of  the  beam-plasma  media.  These 
waves  form  two  orthogonal  sets  in  an  empty  wave  guide,  and  the  approx¬ 
imation  is  that  these  sets  are  still  orthogonal  when  the  wave  guide  is 
not  empty.  To  use  this  approximation,  set 

E0  =  Hr  =  Hz  =  0  ,  (3.32) 

or 

Ez  =  Er  =  Hg  =  0  .  (3.33) 

The  quasi-static  approximation  is  used  when  the  waves  of  in¬ 
terest  are  dependent  on  the  presence  of  a  beam  plasma.  The  approxi¬ 
mation  states  that  the  r-f  magnetic  field  is  small  compared  to  that 
usually  associated  with  an  electromagnetic  waveand,  therefore,  the  r-f 
electric  field  is  nearly  irrotational .  We  can, therefore,  set  the  curl  of 
E  j  equal  to  zero  and  solve  for  the  scalar  potential  <J>  instead  of  using 
a  vector  as  a  variable.  To  use  this  approximation,  set 

V  x  Ej  =  0  ,  (3.34) 

E4  =  -v$  ,  (3.35) 
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and  evaluate  H  ^  from  Equation  (3.  18).  This  is  called  the  quasi-static 
approximation,  because  its  validity  depends  on  the  phase  velocity  of 
the  wave  being  much  less  than  the  speed  of  light.  The  approximation 
is  also  not  valid  at  frequencies  for  which  any  of  the  modes,  e.  g.  ,  the 
beam  cyclotron  modes,  have  a  resonance  resulting  from  the  cyclotron 
motion  of  the  electrons  . 

When  this  approximation  is  used,  the  wave  equation  for  the 
homogeneous  case  {see  Appendix  B)  is 

D  drr  $  +  T  dr  *  +  A  *  =  0  -  (3-36) 


where 


d,  =  — 
1  u> 


pp  pp  pp 

Equation  (3.  3S)  is  a  form  of  Bessel's  equation  and  its  solution  is  given  by 


4> 


(3.39) 


J  (b  ,r  )  = 
o'  oi  o' 


(3.40) 
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(3.41) 


Equation  \3.41)  is  a  dispersion  relation,  whose  solution  gives 


P  -  P(<i,  di.  &2.,  d3) 

where  d,  -  l  A  v  /r  co  .  Figure  4  plots  B/B  versus  a  for  a  set 
3  o  o  pp  a  r  r'  reo 

of  parameters  (d ^  ,  d^  ,  d that  correspond  to  expected  experimental 
values.  Let  us  consider  the  cu*e  for  rQ  approaching  infinity!  then 
bQj  approaches  at-ro  and  Equation  (3.41)  reduce!  tc  the  following  rela¬ 
tions  : 

-P2  u  o  ,  (3.42) 


(3.43) 


(3.44) 
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Equation  (3.  42)  voids  the  quasi-static  approximation;  the  correct  equa¬ 
tion  is 


(3.45) 


and  the  solutions  represent  electromagnetic  waves  that  propagate  in 
free  space.  Equation  (3.43)  gives  the  space-charge  modes, 


Equation  (3.  44)  gives  the  cyclotron  modes, 


a  d1 


The  space  charge  and  cyclotron  modes  are  beam  modes  that 
are  modified  by  the  presence  of  the  plasma.  By  the  physical  argu¬ 
ments  given  earlier,  we  know  that  the  amplifying  wave  must  have  a 
phase  velocity  approximately  that  of  the  beam,  i.e.  ,  (}  Pg  for 

amplification.  This  suggests  the  following  change  of  variable  often 
used  in  microwave  tube  analysis, 


(3.46) 


where  R  is  called  the  plasma-frequency  reduction  factor,  and  "plasma" 
refers  to  the  beam  electrons.  In  our  case,  R  is  usually  complex, 
so  that  the  analogy  does  not  strictly  apply;  nevertheless,  we  will  call 
R  the  beam-plasma  reduction  factor.  Gain  is  proportional  to  the 
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imaginary  part  cf  R  ,  whenever  the  gain  is  possible.  Rewriting  the 
dispersion  relation  in  terms  of  R  gives 


For  r  finite,  the  singularity  at  the  plasma  frequency  in  the  space- 
charge  mode  is  removed.  Indeed  for  the  choice  of  parameters  shown 
in  Figure  4,  neither  the  real  nor  imaginary  parts  of  R  are  large  at 
a  equal  to  cne .  Physically  we  can  think  of  this  as  a  coupling  of  a  grow¬ 
ing  space  charge  mode  to  a.  plasma  wave  that  is  being  attenuated,  with 
the  degree  of  coupling  dependent,  on  the  influence  of  the  cylinder  walls. 
Figures  5,  6,  and  7  show  the  effects  of  gain  versus  frequency  varying 
di  ,  d^  ,  and  d^  respectively.  Whenever  the  gain  is  particularly  high, 
it  is  very  apt.  1,o  be  limited  by  other  effects,  e.g.  ,  inhomogeneities, 
collisions,  nonzero  temperature.  Figure  4  is  the  dispersion  diagram 
of  the  modes  given  by  Equations  (3.  16),  (3.  17),  and  (3.  18).  As  expected, 
the  mode  coupling  through  induced  wall  currents  has  removed  t.he  sin¬ 
gularity  in  the  space-charge  mode  of  the  beam  plasma  near  the  plasma 

frequency  Figure  4  could  have  been  obtained  qualitatively  using  coupled 
19  20 

mode  theory.  ’  The  interaction  of  the  slow  cyclotron  mode  with 

the  forward  mode  leads  to  forward-wave  gain,  and  the  interaction  of 
the  fast  cyclotron  mode  with  the  backward  mode  leads  to  backward-wave 
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gain,  Note  that  by  having  the  inequality  1  > .  ch  >  d  ,  the  two  forward- 

£  i 

wave  regions  can  be  combined  to  get  a  broad  amplification  bandwidth. 
Maintaining  this  inequality  at  high  frequencies,  however,  would  be  im¬ 
practical,  so  this  case  is  not  discussed  further. 


3.  Effects  of  Beam -Plasma  Inhomogeneity 
Using  the  same  restrictions  as  in  the  homogeneous  case,  gives 
differential  equation  (see  Appendix  B): 


Dd 


i=±?l  d  4>  +  A 

y*  -v  1 


4>  =  o 


(3,48) 
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1  +  R 
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(3.49) 
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(3.  50) 
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(3.  51) 


The  density  variations  were  chosen  to  be  of  the  form. 


A.  e 
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The  constants  determine  the  cut-off  points  for  the  variation  while 
the  constants  q^  determine  the  sharpness  of  the  cutoffs.  This  gives 
great  flexibility  in  choosing  density  cross  sections,  as  illustrated  in 
Figure  3.  Examples  of  cross-sections  used  are  shown  in  Figure  8. 
The  coefficients  for  Equation  (3.9)  become 


A  second-order,  ordinary,  linear  differential  equation  with  nonconstant 
complex  coefficients  that  are  also  rational  functions  of  the  eigenvalue  R  , 
is  difficult  to  solve  explicitly  by  analytic  techniques.  It  can  be  handled 
by  a  series  of  straightforward,  though  long,  numerical  computations.  A 
computer  program  for  solving  this  equation  is  given  in  Appendix  C.  The 
program  couples  a  procedure  for  numerical  integration  to  a  search  pro¬ 
cedure  for  R  .  There  was  a  question  of  validity  when  D(r)  =  0  was  in 
the  range  of  integration.  However,  by  expanding  the  coefficients  in  a 
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Indeed, 


power  series  about  such  a  point,  we  obtain  a  good  solution, 
increasing  the  fineness  of  the  integration  at  these  points  gave  finite, 
though  discontinuous,  values  for  <j)  and  d^cf).  Figures  9.  10,  and  11 
show  the  gain  characteristics  for  the  cross  sections  of  Figure  8. 

In  Figures  10  and  11  the  gain  curves  for  the  cross  section  with 

22 

the  smallest  diameter  show  several  peaks,  because  of  lowered  cou¬ 
pling  between  the  beam  and  plasma  modes.  Consider  the  analogy  of 
a  lumped  circuit  of  several  cavities  synchronously  tuned  and  strongly 
coupled  to  each  other  by  a  linear  mechanism.  When  sufficiently  over¬ 
coupled,  the  circuit-frequency  characteristic  shows  just  one  broad  peak. 
As  the  coupling  decreases,  the  peak  becomes  narrower  and  higher  and 
then  becomes  two  broad  peaks.  These  also  become  narrower  and  higher 
until  they  also  divide  into  two  peaks,  etc.  The  results  of  the  repeated 
division  is  shown  in  Figure  10,  where  the  beam-plasma  coupling,  through 
induced  surface  current  on  the  metal  wall  (or  plasma  sheath),  is  reduced 
as  the  diameter  of  beam  and  plasma  is  reduced. 
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FIGURE  1C.  Gain  Dependence  on  Plasma  Width 
(See  Type  b  of  Figure  8.) 


RE  it,.  Cain  Dependence  for  a  HoJ.ov/  Beam. 


into  account,  the  following  analysis  is  accurate  only  to  an  order  of  magni¬ 
tude.  The  solution  for  the  space-charge  mode  in  the  infinite  case  is 


This  removes  the  singularity  at  the  plasma  frequency,  and  if  there  were 
no  other  broadening  processes  available,  e,  g.  ,  finite  geometry,  inhomo¬ 
geneities,  plasma  teinperatur e ,  it  would  give  the  gain  as  a  function  of 
frequency.  Figure  1 1  plots  the  gain  characteristic  fov  different  values 

of  d^  .  Iir  the  plasma  experiment,  we  expect  a  value  for  of  less 

.3 

than  10 


b.  Plasma  Temperature 

Here  we  work  directly  with  the  Boltzmann  equation,  Equation  (3.  21), 

to  express  f^.  in  terms  of  and  v  and  then  integrate  to  get  p  , 

as  shown  in  Appendix  D  and  consider  only  v™  /  0  .  It  is  reasonable 

1  i 

to  consider  only  the  effect  of  plasma  temperature,  since  the  relative 
velocity  spread  of  the  beam  is  much  less  than  that  of  the  plasma.  From 
Appendix  D  we  get 
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FIGURE  12.  Gain  Deperidenc«  on  Plasma  Collisions. 
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The  approximations  are  good  for  all  waves  that  satisfy 


8  «  1 


Consider  that  a  stationary  plasma  (3  .  is  infinite  in  an  infinite  geometry. 


The  dispersion  equation  is 


-?  *♦!(£) 


Figure  13  shows  the  gain  characteristic  as  a  function  of  y,  where 


T  “  v 
o  T 
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C.  CALCULATION  OF  TOWER  GAIN  AND  PHASE  SHIFT 

1.  Experimental 

The  equipment  needed  to  measure  power  gain  is  shown  in  Figure  14a 


A  reading  was  taken  on  the  calibrated  attenuator  with  the  amplifier  removed 
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(dashed  line).  The  amplifier  was  then  placed  in  the  circuit  and  the 
attenuator  readjusted  to  get  the  same  reading  on  the  power  indicator. 
The  difference  in  attenuator  settings  represented  the  power  gain.  Of 
course,  corrections  resulting  from  mismatch  or  amplifier  accessories 
had  to  be  taken  into  account. 

In  Figure  14b,  equipment  was  set  up  to  measure  amplifier  phase 
shift  to  within  an  additive  constant  of  2k-rr  radians.  A  minimum  was 
found  on  the  slotted  line  with  the  amplifier  out  of  the  circuit  (dashed 
line).  With  the  amplifier  in  the  circuit,  a  second  minimum  was  found  . 
The  amplifier  phase  shift  was  equal  to  twice  the  VSWR  minimum  shift. 


2.  Theoretical 

To  compare  the  theoretical  and  experimental  values  of  power 

gain,  we  must  know  the  coupling  coefficients  for  the  input  and  output 

* 

klystron  gaps.  The  procedure  for  determining  these  coefficients  is 
lengthy,  and  we  shall  assume  that  they  are  known.  The  total  power  gain 
is  given  by 


E2 

3ei 


2 

**11 
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(3.  59) 


where  p^  is  the  ratio  of  the  power  of  the  lowest-order  mode  leaving 
the  input  cavity  to  the  input  power;  p^  is  the  ratio  of  the  power-of  the 
lowest-order  mode  in  the  output  section,  to  the  power  output  from  the 
amplifier;  and  Vg(R)  is  the  velocity  of  energy  transport,  which1  we  take 
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equal  to  the  group  velocity  of  the  wave.  For  the  cases  discussed  in 
Ill  B, 
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III.  CCfNCLUSIONS 


We  have  found  that  the  effects  of  finite  geometry,  inhomogeneous 
density  cross  sections,  plasma  collisions,  and  plasma  temperature  all 
limit  the  gain  and  broaden  its  frequency  range.  For  the  experiment  now 
being  conducted,  we  expect  the  density  cross  section  to  be  the  most 
important  of  these  factors.  For  higher-frequency  applications,  the 


effects  of  plasma  collisions  may  dominate  the  broadening  processes 

1  /2 

since  vc/u  varies  as  n  '  This  is  a  strong  argument  for  plasma 
generation  using  contact  ionization,  which  has  a  percentage  of  ionization 
greater  than  50  percent,  for  frequencies  over  100  Gc/s  . 
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APPENDIX  A.  BOLTZMANN  AND 
BOLTZMANN  TRANSPORT  EQUATIONS 


1.  Boltzmann's  Equation 

Consider  the  distribution  function  f.  .  This  specifies  the  num¬ 
ber  of  particles  of  type  i  at  a  given  point  (x.  ,  v)  in  phase  space. 

U  sing  Liouville's  theorem  that  the  total  de rivative  of  the  distribution 
function  along  a  flow  line  in  the  phase  space  is  zero,  we  get 

f.(t)  =  f.(t  +  dt)  (A.  1) 

Including  the  effects  of  collisions,  we  get 

/Of.  1 

f.(x  1  *  tU  •  ^  +  v  dt  ,  t  +  dt)  -  f .  (x  ,  v  ,  t)  =  I  dt  (A.  2) 


Expanding  f^  (x  ,  v.  >  1)  about  t  ,  gives 
Sf.  9f .  Ox.  Of.  Ov. 

f .  (t)  +  ■  dt  +  g—  •  -g-j—  dt  I  g v'  ’  -g—  dt  +  higher  order  terms  in  dt 

—  i  —  i 


(A.  3) 


Collecting  terms  of  order  dl  and  letting  dt  "*■  0  .gives  Boltzmann's 
equation , 


V  +  a 


(A.  4) 


-50- 


2.  Boltzmann's  Transport  Equation 


When  the  time  scale  is  greater  than  the  time  between  colli¬ 
sions,  the  particle  density,  average  velocity,  and  energy  can  be 
treated  directly,  i.c.  ,  without  velocity  distribution  functions,  since 
these  quantities  are  invariant  during  a  collision.  Integrating  Equation 
(A.  4)  over  velocity  space,  we  get 

ur  =  -v  ■  +  \  ■  (A-v 

where  L  represents  net  density  of  particle  creation,  and  v^  is  an 
average  velocity  defined  by 


v 


1 


f .  (u  .)  d" 

i—i 


(A.  6) 


Tlie  right-hand  side  of  Equation  (A  5)  represents  the  net  particle 

flux  from  a  differential  volume  element.  Contributions  to  this  term  can 

arise  in  three  ways:  (1)  as  a  result  of  diffusion  in  co-ordinate  space-, 

+  V  '  5  n  ’  (2)  as  a  result  of  particle  mobility  in  an  electric  field 

(akin  to  a  diffusion  type  process),  -  V  ■  p.  •  E  n  ;  (3)  as  a  result  of  dif - 

djE  v  n 

fusion  in  velocity  space,  -  -  and  -  •=—  •  (Note  that  these  two 

m  v  L 

terms  are  usually  small  compared  to  the  rest,  unless  the  mean-free 
path  is  greater  than  or  approximately  equal  to  a  typical  dimension  of 
the  geometry). 
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APPENDIX  B.  DERIVATION 
OF  POISSON'S  EQUATION 

Using  the  quasi-static  approximation,  =  0,  the  electric 

field  can  be  derived  from  a  scalar  potential.  The  set  of  equations  to 
be  solved  is: 


E  i  =  -  V< 


E  ,  =  p,  /  e 
—  1  vo 


dtnu  +  V  •  n0llil+  nn^01  =  0  • 


Vl2+v‘  (\uXl2  +  n12X02/ 

Vll  +  (l0l'V)  —11  +  (in 


X01 


i (1 1  +  In  x  2.0  +  1  o i  x  —  l 


dtll2  + 


(xoz  • v)  Xu  +  (Xu-v) 


-  r,  E1  +  V12  x  Bo  1  V0(,  x  P4 


(nll  +  nu) 
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The  following  conditions  apply  to  the  beam-plasma  problem: 


»1 

=  0 

V  '  —02 

=  0 

dt 

=  jw 

Xoi 

=  0 

E0 

-  0 

d 

z 

=  -JP 

X02 

A 

=  kV0 

de 

=  0 

T> 

±±0 

1! 

w 

o 

Using  Equation  (  B.  1)  in  Equation  (B.  2)  gives 


(V  •  V)  <|> 


(B.9) 


The  y  in  cylindrical  co-ordinates  is  written 

V  ^  =  d4,^+~-dQiV0+<iiVk 

-r  r  ~  r  6  .2 

(V  -V)^  ~  +  4  d0O  +  +  + 

r 

v  •  A  —  d  (rA  )  l  —  d  A  t  d  A 
—  r  r  '  r'  0  0  z  z 

U  sing  Equations  (  H .  H) ,  (11.10)  in  (13.1),  (11.4),  ( 11 .  5) ,  (11.6).  and 

(  U.  9),  we  get 

jwiijj  1  v1irclrnoi  +  noirdr  (r  Vlli-)  1  “ol'"J^V  1  lz  0  • 

■’w"t2  1  V  12rdrn02  +  n027  <MrV12r)  +  no  1^ V  1  2z  +  V0 '^n12 

JW X  i  i  =  -  n  +  Vu  x  B0j 

12  +  Vo  —  12  =  -  n  v  ^  +  X  12  x  — o)  ’ 


(B.  10) 


(B.  1  1) 

0  , 

(  13.  12) 
(B.  13) 

(  B.  14) 


-hi- 


(B.  15) 


Next,  we  solve  for  the  small-signal  velocities  in  terms  of  the  scalar 


potential.  Expanding  Equation  ( B.  14),  we  get 


(B.  16) 
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The  following  identifications  are  made: 


“c  =  ^B0 


U) 

v7 


D  =  1  - 


1  - 


Following  a  similar  procedure  for  Equation  (  B.  13),  we  get 


■  rj 


0 


d  <b 
r  T 


-Jco 


Wc 
1  _  -it 


-jP* 


•  11 


jw  1  - 


(B.  17) 


Substituting  Equation  (  B.  16)  and  (  B.  17)  into  Equations  (B.  11)  and  (B  .  1Z), 
we  solve  for  the  small-signal  electron  densities.  Expanding  Equation 
(B  .  1 Z) ,  we  get 


/.  \  (dr4>)  (drnn?) 

jwr)12  +  (js?)  - 7T -  +  Vq  (-jP)  n  2 

'  K )° 


+  p2 

\}u 


(-j^)(-jP)  nQZ  /  >  np 2.  (drr*7dr)  * 


i1' 


JL 

Pe 


(zn\ 

v10  / 


=  o  , 


1  -  -it-  D 
Pe/ 
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12 


(B.  18) 


n 


P  “0  2* 


(drn02){dr<t») 

n02  (drr  +  7dr)  + 

D 

D 

(ju)* 


2 


Following  a  similar  procedure  with  Equation  (B.  11),  we  get 


‘11 


n 

(jw)Z 


P2n 


01 


4>  - 


(drn01)(dr<W 


(B.  19) 


Substituting  Equation  (B.7),  (B.  18),  and  (B.  19)  into  Equation  (B.  15), 
we  get 


-  f?6- 


-  r>7  - 


(a)  , 

J£b 

U) 

PP 

A 

n0i 

Hpi 

noi 

max 

V 

o 

A 

*>02 

^02 

r  w 
o  pp 

^02 

max 

(B.  21) 
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APPENDIX  C.  COMPUTER  PROGRAM 
FOR  THE  INHOMOGENEOUS  CASE 


i.  Derivation  of  Parameters 


The  differential  equation  is 

"D(R  .  v)  +  E(R  v)l 

D(R,  v)  dvv4>(v)  +  I - - - - - i-jdv<t>(v)  +  A(RlV)+(v) 

subject  to  the  initial  conditions, 

<M°)  =  0 

dv<t>(o)  =  0 

The  complex  parameter  R  is  to  be  chosen  such  that 

<j>(v)|  =  0 

|v=l 

Expanding  in  real  and  imaginary  parts  we  get 

4>  =  +  j  <t>2  , 

D  =  D1  +  jDa  , 

E  =  Et  +  jE2  , 

A  -  A^+jA^  , 

R  =  Ri  +  jR2 


0  , 

(C.  1) 


(C.2) 


(C.  3) 


Substituting  Equation  (  C.  3)  into  Equation  ( C.  1),  we  get  a  system  of 
two  real,  coupled,  linear,  second-order  differential  equations: 


-59- 


t  F  (R^  «  .  <*v4>j(v)  I  Rz,  v)  dv<(>2(v) 

•)  >R ,  ( v )  I  FS(R},  R2.  v)  4>2(v)  0  ,  (C.4) 

t  (  I  .  r  K2,  vJd^M  -  F  ^(R  j,  R2,  v)dvi))1(v) 

'■>  •:*><*■)  '  Fs(.<1.H2,v)4,1(v)  =  0  .  (C .  5) 

1  ii  i •. u u»  <tie. 

10  <|>2(o)  0 

(i  il  io  ,(u)  0 

a  i)  o  .  (c.  6) 


<VoV) 


D1  '  1  - . -<T 


fR  ,  p 

1 


'-I 


D2  =  CR?_  e 


2  j2 

C1  =  Q  -  dl 


+  CR2Ae 


•  ,r  v) 
3  u 


R?  -  Ri  ■  4 


*5  -  P4  -  *  (w^zY 


{<•  HI) 


(;,<v4ai  rov> 


('l  <;VnV) 


1  1  5A  (‘'  ')  n"' 


H,  <" 


T'V.  CR4(;'Rr<»V)  " 


‘1?.  ■•<VuV) 


I  .  ■  ■) 


1  •'  ' 


l,  I 


ql 

CR33 

a2  -  d2 

cr3  = 

«1 2CR1  1 

CR3A  = 

<13CR1A 

1 

(C.  13) 

-(a.  r  v)  1  -(a?r  v)  ‘(a3rr>v) 

A,  -  CRS-CR„,  -CR567e  '  CR567Ae 


{ C.  14) 


A,  =  CR7-CR73e 


-{a.r  v) 
'  1  o 


CR765e 


-(art) 
C  o 


CR765Ae 


(a3rot) 


(C.  15) 


CR567A  PCR567 


-f,Z- 


CR^R^)  -  CR^Rj  -  R2) 


!R1  -  Rtf  4  K^f 


CR765A  =  PCR765 


(C.  16) 


The  integration  procedure  requires  that  the  real  and  imaginary  parts 
of  dvv.(jj(0)  be  known.  To  find  theBe  we  write  a  power  series  expan¬ 
sion  in  v  about  the  origin: 

<|>{v)  =  1  +  (a^  +  ja^)  t2  +  ...  + 


(v)  =  2(at  +  ja^t  +  ...  + 


Mv)  =  2(at  +  ja2)  +  ...  + 


(C. .  17) 


dvv<t>(°)  =  2a  d  +  jZa2 

To  find  a^,a2  we  substitute  Equation  (C.  17)  into  Equation  (C.  1): 
2(2ai  +  2ja2)  -  CRj  -  CRJA  +  j(CR2  +  CrJ 

1  J 

+  [cr5-cr53-cr567  -cr567A  +  (CR7-CR73 


'CR765  "CR765A5] 


[c°N1  (i-JL  -CRj-CR^j  *  CON  2  (cR^CR^j) 


4  1  -  K  -  CR1  *  <CR2  *  CR2A|Z 


(C.  18) 
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CON  2 


M-cr‘) 


CONI  CR2  +CR2A 


i  \  l  \ 

--  -  CR, ')  t  CR,  t  CR,  ' 

V  \  2  2AJ 


con  i  =  -cr5+ck53  +  cr567  +  cr567A 


CON  2  =  =CR7+CR73  CR765  -CR765A 


(C.  19) 


2.  Computer  Program 


a.  Elements  of  Program 

The  computer  program  consists  of  seven  subprograms, 
t.  UELMAt; 

Tl-is  is  the  main  program.  It  reads  the  input  data  and  com¬ 
putes  ()u*  eniiKi a nt.s  that  are  independent,  of  R.  This  program  also 
regulates  the  iterative  search  procedure  for  the  value  of  R  that  min- 


iniiy.es  |<|>(1)|  .  Starting  from  a  given  value  of  R,  it  computes  in¬ 

crement  s  Ae.  until  |  <{>(  1 )  j  is  minimized. 


/.  RUNC1E 

This  is  a  standard  routine  designed  to  integrate  a  system  of 
differential  equations  by  the  Rurige  -Kutta  method.  First  Lhe  system 
is  rerlneed  In  Lhe  following  form: 

rlv  y  ^  ( v )  f^Hj  ,  7.? . v) 


,lv'/n(v)  fn(x'l  ’  */' 


V  v) 


(C.  2d) 


For  on  i-  case  l)  1-  K^na'imi  (i  . '))  hm  umes 


d  7.  ~' 

d  A, 

v  l 

V  1 

d  Z  -v  = 

V  2 

d 

v  rZ 

dvZ3  " 

d  <h 

vv  1 

d  z  .  = 

d  6™ 

v  4 

vv1  2 

-  jV^z^  +  F^z^ +  F^z^ +F^  •  z^J  /F^  , 

"  [F2Z4  "  F  3Z3  "  F4Z  2  "  F  52 1]  /  *’  1 


(C.  21) 


Given  the  starting  values  Z^(0)  ,  this  routine  integrates  Equa¬ 
tion  ( C.  ?.\)  from  v  -  0  to  v  ~  1.0  .  During  a  typical  step  in  this 
program,  say  from  t  to  t  I  A  ,  the  routine  branches  out  to  another 
subprogram  called  DER1V,  which  computes  d^z^  ,  four  times.  During 
these  limes,  the  values  of  z.  are  stored  in  temporary  storage. 


1.  J5F.R1V 


This  ituhprugram  on  1  cn  I  at  e  s  the  derivatives  d  z.  an  given  in 

Equation  (G.  21) .  If  v  0,  then  f  ^  is  set  equal  to  d  v<j>^(o)  2a  ^ 

and  f,  is  set.  equal  to  d  Ado)  since  in  the.  normal  evalua- 

4  vv  ■  f.  A 

lion  of  d  it,  there  in  a  division  tv  v  . 
v  v  1  ' 

A  test  in  made  on  the  variable  N  (the.  number  of  equations 
in  the  system),  whii  )■  is  usually  equn  I  to  four,  if  the  magnitude  of 
i'\j  ii  less  t!i.‘n  !0  '  ,  Hie  value  of  N  v/ ill  he  set.  equal  to  two.  For 


H  equal  to  t v/* > ,  a  reduced  system  of  equations  is  integrated  (see 


4.  EXPRESS 


This  subprogram  is  entered  from  RUNGE  when  v  -  0.0  and 
=  1.0  .  When  the  value  for  R  that  minimizes  |<|>(l)j  is  found, 
this  program  is  entered  for  all  v  ,  so  that  every  step  in  the  integration 
is  printed  for  this  final  calculation. 

The  main  function  of  thiB  program  is  to  print  a  line  of  output 
values  of  v^  ,  <t>2,  d^,  dv<j>2,  d^^,  eivv4»2’  Dv  D2  at  the  proper 
values  of  v  .  It  also  tests  the  magnitude  of  F^  against  10  .  If  it 

is  less,  then  Equation  (C.  21)  is  reduced  to  the  following  svstem  of  two 
first-order  equations; 


dv+l  =  -  +  P2(i52)/,p3  1 

dv^2  =  -(Pl<l>2-  P2^P3 
P1  "  F2F4  +  F3F5  > 

P2  -  F2F5  -  F3F4  ' 

P3  ~  F2  +  F3  • 


(C. 22) 


(C.  23) 


This  system  is  integrated  over  the  next  interval  Av  .  The  interval 

is  divided  into  ten  steps  of  length  Av/10  ,  and  after  every  10  steps 

the  magnitude  of  F^  is  tested  again.  When  the  magnitude  of  F^  is 
~S 

greater  than  10  ,  then  N  is  reset  to  four  and  the  integration  inter¬ 

val  is  reset  to  Av  . 

6.  PHICON 

This  routine  is  entered  from  the  main  program  DEIMAG.  Its 
function  is  to  set  up  aLl  the  necessary  starting  conditions,  enter  RUNGE 


-66- 


to  integrate  Equation  (C.  21)  for  a  particular  value  of  R,  and  to  com¬ 
pute  <b(l)^  .  It  also  tests  |tj>(l)|  to  see  whether  the  search  procedure 
should  be  stopped. 

6.  RCON 

This  is  an  auxiliary  routine  entered  from  PHICON.  Given  a 
particular  value  of  R  ,  it  computes  all  of  the  R -dependent  constants 
used  in  computing  the  derivatives  (namely  CR^,  CR2’  CRo>  .  ■  •  ,  etc.  ) . 

7.  TCAL.C 

This  is  an  auxiliary  routine  entered  from  DERIV.  For  a  par¬ 
ticular  value  of  v  ,  it  uses  the  stored  R-dependent  constants  and  com¬ 
putes  D^,  D^,  E ^ ,  E-,,  A^,  A^,  F ^ ,  F^,  F^,  F^,  F^. 

b.  Search  Procedure 

An  initial  choice  for  R  is  given  the  computer  along  with  all 
other  input  data.  From  this,  four  alternate  starting  values  for  R  are 
calculated  and  stored  for  future  reference.  Get  A^  be  the  magnitude 
of  <|)(1)  squared  after  the  i^'1  integration.  PHICON  supplies  A.  to 
DE1MAG,  but  first  tests  it  against  .01  .  If  it  is  smaller,  DE1MAG 
directs  EXPRESS  to  print  the  output;  if  it  is  larger,  DE1MAG  stores 
A.,  changes  R  by  A,  compares  A.  +  ^  with  A^  .  If  A^  (  j  is  larger 
than  .01  ,  the  increment  in  R  changes  direction  by  tt  or  ir/2  radian 
gets  larger  or  smaller  dapending  on  the  previous  values  for  A  and  A 
The  search  ends  when  either  A  is  so  small  that  increased  accuracy 
would  not  justify  the  added  computer  time, or  when  A^  is  less  than  .01. 


Superimposed  on  the  above  procedure  is  an  additional  routine 
that  prevents  the  computer  from  searching  a  nonconvergent  region  of 
the  R  plane.  In  working  with  the  program,  it  was  noticed  that  if 
l^v)!  was  greater  than  one  near  v  =  o,  no  good  values  of  R  were 
found  in  that  region.  So,  a  test  is  made  during  the  integration  pro¬ 
cedure.  If  j { •  0 Z) j  is  greater  than  one  for  two  consecutive  values 
of  R  ,  DE1MAG  selects  an  alternate  starting  value  for  R  that  was 
held  in  storage  from  an  earlier  part  of  the  program.  DEIMAG  will 
go  through  this  alternate  procedure  four  times,  if  necessary.  The 
fifth  time,  it  will  go  on  to  the  next  input  data  card.  The  flow  chart 
shown  in  Figure  C.  1  for  DEIMAG  gives  a  more  sequential  explanation 
of  the  search  procedure. 

c.  Input  Preparation 

There  arc  two  programs  using  this  Equation  (C  .  1).  The  first 
is  called  the  "regular"  program,  and  is  the  one  with  the  automatic 
search  for  R.  The  second  is  called  the  "grid"  program,  and  does 
not  have  the  automatic,  search  for  R  .  Instead  it  takes  a  rectangular 
grid  in  the  complex  U  plane  and  transforms  it  into  a  corresponding 
area  in  the  complex  <|>  plane. 

The  regular  program  has  two  types  of  input  data  cards,  Card 
types  R-l  and  R-2. 

Card  Type  R-l 

Columns  1-7  8-14  15-21  22-28  29-35 

Data  d,  d,  d„  a,r  i,r 

1  2  3  1  o  2  o 
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Columns 


36-4?.  43-49  50-56  57-62 


63-69 


Data 


q?  ^3 


P 


C  a  r d  Type  R  -  2 


Columns  1-10  11-20 

Data  a  R, 

1 


21-30  34-35 

R2  1,  -1,  0 


The  values  in  columns  34  and  35  of  card  type  R-2  indicate  pro¬ 
cedure.  No.  1  means  that  the  next  input  data  card  must  be  of  card 
type  R-l.  No.  0  means  that  the  value  at  R^,  R^  on  this  card  must  be 
used  as  the  starting  point  for  R  .  No.  -1  means  that  the  values  of 
R  ^ ,  R ^  derived  from  the  last  input  data  card  will  be  used  as  the  start 
lug  values  for  R_j,  R^  .  The  grid  program  also  has  two  types  of  input 
1 1- 1 1 .  i  iiir'ls,  the  first  of  which  is  identical  with  i  a  id  type.  R-l;  the 

si  i  uinl  is  r  .1  id  type  (  I 
fird  Type  ( i  2 

‘  lo I lue.ic.  1  II)  11  20  21-10  il  40  11  50  51  60 

P  it. i  |.  I<  (  If  DCs  I.X1Y  IT'S 

1  el  lyjie  G  ?  “il!  j  m  ■  i  u  i  .ile  n  grid  in  I  In-  K  plane  whose  corners 
a  *  upy  the  l.dlowing  eo -o  rdi  nates  i 

(L<  ,  K  >  ;  (R1,  H  ^  +  1)GY  X  EPS)  ;  (K^  1  DGX  x  EPS);  (R^l  DGX  x  EPS, 
R2  l  DGY  x  EPS)  . 


00- 


PROGRAM  DEI  MAC. 

STORE!  *  l 00.0 
ST0RE2  =  1 00 • 0 

DELTA  »  0.01 

1  READ  lOO.Sl .S2.D3.AlR0.A2R0.Ol .02 
100  FORMAT  17F10.6) 

6000  READ  6001 .SETA .A3R0. 03 

6001  FORMAT  OF  10.6) 

999  FORMAT  (13) 

PRINT  3000*  SI  .S2.D3. A1»0, A2R0.Q1  ,02 
3000  FORMAT (IX* 7F 1 0. S ) 

SETA  =  0.0  -  SETA 


PRINT 

605 • BET A  «  A 3R0 • 03 

605 

FORMAT 

(2X.3F10.6  1 

IF  (SI 

-100.0)  555*556,566 

556 

STOP 

555 

CONT I NUE 

FL  -  1.0/DELTA 
2  READ  1 01 . ALPHA ,R1 *R2* I M 
KC  =0 

101  FORMAT (3F 10.5. 15) 

PRINT  3C01.  ALPHA 
3001  FORMAT  <  8H  ALPHA  =.F10.5> 

calculations  for  r- independent  CONSTANTS 


IF  (III) 

3005*2003, 1 

3005 

R 1  =  STORE  1 

R2  =  ST0RE2 

STORE t  * 

l  00.0 

S  T  ORE 2 

a  100.0 

2005 

IF  (Rl  - 

100.0)  2003,2.2 

2004 

STOP 

2001 

CONTINUE 

RAL)  a  Rl*f?l  4  R2»R2 
ADD  =  0.01*  SQRTF  (RAL) ) 

PF  ST  «  0.01  »ADf) 

RISE  r  (  I  I  *  R 1 /2 • O 
R l SET ( 2  )  *  R 1 »?  .0 
R2SET(  l  )  -  R2/2.0 
R2SFT ( 2  )  ■  R2*2.0 
PRINT  900.  R2SCT ( P ) 

990  FORMAT  C  F 1 0.5/  ) 

LBACK  =  o 
I  SET  =  I 
Ho  1.0 

333  RR ( 1 )  -  R I 

RR(?l  =  H*» 

WAU  :  ttl*R|  ♦  R2»W2 
RAO  =  SQRTFCRAH) 

DHL  I  “  p  AO*  0.01 
DHL 2  «  DELI 
KC  «  KC  4  1 

CAUL  PH  ICON 
IF  ( LBACK 1  21  *  301  .300 

301  CONTINUE 

COMP  e  PAD 
KEY  «  1 

105  DL  <11  -  DEL  l 

OL  (  2  >  *  DE  L  2 

1  0  1  RR(KEY)  -  RR  (KC  Y  )  ♦  OL  (KEY) 

R  1  a  RP ( 1  ) 


R 2  =  RR ( 2 ) 

KC  =  KC  +  1 
PRINT  999, KC 
IF  (KC  -  42)  307.300,300 

307  CONTINUE 

CALL  PH ICON 
IF  (LBACK)  21*303,300 
303  CONTINUE 

!  F  (PAO-CO^P)  14,11*11 

11  RR(KEY)  *  RRCKEY)  -  2.0*DL(KEYI 

DL (KEY  >  =  -DL(KBY) 

R  1  8  OR ( 1  ) 

R2  =  RR ( 2 ) 

DELI  =  DL ( 1 ) 

DEL2  =  OL ( P ) 

KC  °  KC  *  1 
PRINT  999. KC 
t f-  (KC  -  42)  306,300.300 

306  CONTINUE 

CALL  PH  I  CON 
IF  (LBACK)  21,304,300 
304  CONTINUE 

|F  (PAD-COMP)  14,14.13 

13  RR(KEY)  *  RR(KLY)  -  DL<KEY> 

GO  TO  200 

14  TEST  =  ( COMP— PAL)  > /DL  (  Kt.  Y  ) 

COMP  =  PAM 

IF  (TF5T-H)  100*171*101 

100  TDEL  =  2  •  0*  AHSF  (  DL  ( K(-  Y  )  1 

IF  (  TDEL  — RAD  )  10T.lO?.IO;’ 

103  DELI  =  2 • O*  DE  L 1 

DEL 2  =  2.0#DIL? 

GO  TO  105 

10?  DLL  1  n  (1.01  *RAR 

DFL2  r  DM  l 

10  IF  (KEY- 1)  106.107*106 

106  KEY  =  1 

GO  TO  106 

1 0  7  KEY  a  ? 

GO  TO  1  PM 

20 0  DL  (KEY  )  =  DL  (Kl- Y  )«  ( -0  .  1) 

T  ST  1  »  AHSF  <D1  (  1)  ) 

T  ST  2  =  AljSF  (PL  (2)  ) 

DHL  1  »  DL (  1  ) 

DLL  2  *  DL (2) 

|F  ( MAX  IF  ( T6T  I  . T S  12  )  -  PL 5 T  )  2 1  . 2 l  ,  1 O 

21  L  ■  I 

KC  =  KC  4  1 

CALL  PHICON 

305  STORF 1  -  R 1 

SfORl-2  =  R2 
GO  TO  2 

300  R 1  =  Rl SET (  I  SET ) 

R2  -  R2SET  <  I  SET ) 

KC  »  O 

] SET  -  !SFT  +  l 

IE  (| SET -31  302,302*2 

3C2  GO  to  J33 

O  STOP 
END 
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APPENDIX  D.  CORRECTION  TERM  FOR 
MAXWELLIAN  DISTRIBUTION  IN  LONGITUDINAL  VELOCITY 


We  can  treat  the  effects  of  nonsingular  longitudinal-  and 
transvers e -velocity  distributions  separately,  since  these  modes  are 
not  coupled  to  each  other  in  Boltzmann's  equation.  From  Equation  (21), 
evaluating  the  contribution  of  the  axial  variation  of  velocity  to  the  par¬ 
ticle  density  variation,  we  get 


oo 

nl  -  /  fl 


(vz)dvx 


- — —  ~  If  (v  )]  dv  .  ( D.  1) 

,(td  -  jlv^  9vz  l  °  z 


For  a  Maxwellian  distribution, 


f  (v  ) 
o'  z' 


I  ,  V 2 
-i/2  -i  ’  Vz"vo  /VT 
it  v  e  / 


,  e  .  -1/2  _  -3  f 

n.  =  (-2  — jtt  <  E.  v_ 

1  m  J  lz  T  J 


.  ,2/  2 
CO  ,  -  Vz"Vo  /  VT 

r  (vz-vo)e  / 


(u  -  (3vz) 


where  v„  =  (2kT/m)  '  .  If  we  set 


F  =  jfl  ir 
P  Jl 


Vz  -fVz-Vo7VT 
-  e  /  d  v 


<n-pvz 


I  mP 


[JvT  "Wo'] 
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Stix  '  has  evaluated  expressions  of  this  form  using  conloih 
integration.  Using  his  result, 


( n.  4) 
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PLASMA  DIAGNOSTICS 


D .  M .  Kerr 


ABSTRACT 


This  report  describes  the  design  and  use  of  singly  re-entrant 
cavities  in  an  electron-beam  hydrogen-plasma  system.  The  cavities 
serve  as  the  input  and  output  coupling  structures  of  a  beam-plasma 
amplifier  and,  in  addition,  they  are  used  to  investigate  the  properties 
of  the  plasma  discharge.  Perturbation  of  a  resonant  cavity  by  a  plas¬ 
ma  is  discussed  and  equations  arc  obtained  relating  the  shifts  in 
resonant  frequency  and  cavity  Q  to  the  average  electron  density  and 
collision  frequency.  Some  typical  experimental  results  are  given  and 
a  linear  relationship  is  found  between  resonant  frequency  shift  and 
discharge  current. 
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1.  INTRODUCTION 


Two  singly  re-entrant  cavities  of  the  type  ordinarily  used  in  kly¬ 
stron  amplifiers  were  designed  and  tested  for  use  in  the  study  of  an 
electron-beam  hydrogen -plasma  system.  Their  primary  purpose  was 
to  serve  as  the  input  and  output  coupling  structures  for  a  beam-plasma 
amplifier.  It  was  also  thought  that  they  might  be  used  to  diagnose  the 
properties  of  the  plasma  and,  in  particular,  to  determine  the  electron 
density  and  collision  frequency.  The  following  sections  describe  the 
design  and  cold  testing  of  the  cavities  and  their  use  as  diagnostic  tools; 
and  presents  the  experimental  results  obtained. 


II.  CAVITY  DKSIGN 

The  beam-plasma  system  operates  at  about  3000  Mc/s  with  a 
5  -kv  electron  beam.  These  specifications  together  with  certain  struc¬ 
tural  details  of  the  vacuum  system,  fixed  the  cavity  dimensions.  Fig¬ 
ure  1  shows  the  quantities  considered.  The  inside  diameter  of  the  drift 
tube  r  was  fixed  by  the  si/.e  of  the  hollow  cathodes  used  to  generate 
the  plasma,  which  were  0.500  in.  in  diameter.  The  stainless  steel  rods 
that  supported  the  plasma  cathodes  and  cavities  placed  a  restriction  on 
the  outside  diameter  of  the  cavities  of  2.50  in.  The  copper  tubing  tor 
the  outer  cavity  wall  had  an  outside  diameter  of  2.00  in.  and  an  inside 
diameter  of  1 . 75  in.  ;  that  for  the  drift  tube  bad  an  outside  diameter 


FIGURE  1.  Cavity  Co-ordinates. 

of  .  750  in.  and  an  inside  diameter  of  .625  in.  Transit  angle  and 
resonant  frequency  requirements  were  then  used  to  find  the  appro¬ 
priate  values  for  h  and  d  . 

The  transit  angle  is  given  by 


-2- 


where  f  is  the  resonant  frequency,  and  where  the  beam  velocity  is 


u 

o 


It  was  desired  to  keep  8^.  between  1  radian  and  2  radians  in  order  to 
assure  a  reasonably  strong  electric  field  across  the  gap.  This  meant 
that  d  had  to  be  between  0. 0876  in.  (0^,  =  1)  and  0.  1752  in.  (0rj,  =  2)  . 

A  second  consideration  in  establishing  the  gap  length  was  the  problem 
of  multipactor ,  or  breakdown,  in  the  gap  caused  by  secondary  emission 
in  phase  with  the  gap  voltage.  A  first-order  treatment  of  this  problem 
showed  that  it  would  not  be  a  factor  at  the  power  levels  contemplated 
for  the  system.  The  problem  was  approached  by  considering  an  elec¬ 
tron  with  7.ero  velocity  at  one  side  of  the  gap  when  a  r-f  voltage,  is 
applied.  The  force  on  the  election  is  given  by 

V 

e  K  =  e  — sinwl  =  m  a(t)  ,  (2) 


and  since  the  acceleration  is  known,  the  position  can  bo  found  by  inte¬ 
grating  twice  and  using  the  initial  conditions  u(l  c  0)  =  0  and 
7,(1  :  0)  -  0  ,  which  gives 
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For  multipactor  to  occur,  the  electron  must  arrive  at  the  other  side 
of  the  gap  (that  is,  z  =  d)  when  tot  =  tt  (i.e.  ,  just  as  the  phase  of  the 
r-f  voltage  reverses)  with  sufficient  energy  to  cause  secondary  emis¬ 
sion:  i 


e  V 

z  (<ot  =  ir)  =  — = — —  it  .  (5) 

<o  md 


To  get  an  idea  of  the  level  at  which  multipactor  effects  could  begin  to 
appear,  one  can  solve  Equation  (5)  for  V_^  : 


I  <- 

-  ^L.  — 
[  TT  “ 


(6) 


Solving  Equation  (6)  for  several  gap  spacings  gives  the  following 
results: 


d  (in.  ) 

0^,  (radians) 

V  (volt  a) 
m 

0. 0876 

1. 00 

3190 

0.  100 

1.14 

4160 

0.  125 

1 . 43 

6500 

The  maximum  value  of  r-f  voltage  across  the  gap  is  V  ,  and  the 
experiments  planned  for  the  beam-plasma  system  should  not  require 
values  of  V  greater  than  those  listed. 

The  next  step  in  the  design  procedure  was  to  consult  the  design 
charts  in  Moreno*  to  determine  the  values  of  h  and  d  that  would  re¬ 
sult  in  a  resonant  frequency  of  about  3000  Mc/s  .  Several  possible 
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pairs  of  h  and  d  values  were  chosen  and  were  then  checked  using 
the  design  equations  of  Fujisawa.  There  was  about  a  10  per  cent  dif¬ 
ference  between  their  predicted  resonant  frequencies.  After  the 
cavities  were  cold  tested.  Fujisawa's  equations  were  found  to  give  the 
best  agreement.  The  sample  calculation  using  the  dimensions  of  the 
cavities  built  for  ro!d  testing  follows: 

r  ^  =  0  875  in.  =  2.  22  cm  d  =  0 .  1  25  in.  =  0.318cm 

r  ^  =  0.375  in.  =  0.952  era  h  =  0.500  in.  =  1.280  cm 

r  -  0.3125  in.  =  0.795  cm 
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Zu 


1 


3260  Mc/s 


(8) 


°  “  n/LC 


■J  MCj  +  co) 


Two  cavities  were  built  for  cold  testing;  resonant  frequency, 

Q  ,  and  R  /Q  were  determined.  Tlie  test  methods  are  discussed  in 
s  o 

the  following  sections  since  the  calculations  are  important  in  the  Later 
use  of  this  type  of  cavity  for  plasma  diagnosis.  The  cold-test,  results 
led  to  a  few  minor  modifications  in  the  design.  The  cavities  to  be  used 
m  the  actual  system  were  then  built  from  the  design  shown  in  Figure 
2.  Figure  3  shows  the  method  of  loop  coupling  and  the  vacuum  seal 
in  the  coaxial  input  line. 


III.  MEASUREMENT  OF  CAVITY  Q 

Cavity  Q  was  measured  using  the  standard  point -by -point 
impedance  technique  as  described  in  Gnixton.^  When  calculating  the 
Q  of  resonant  cavities  for  use  in  plasma  measurements,  the  neces¬ 
sary  accuracy  is  not  achieved  if  one  neglects  the  series  losses  in  the 
coupling  circuit  between  the  cavity  and  transmission  line.  Since  the 
total  loading  of  the  cavity  determines  the  width  of  the  resonance  curve, 
the  loaded  Q  best  represents  the  effective  load.  For  diagnostic  work 
the  various  Q  values  are  computed  so  as  to  include  the  effects  of 
colliding  losses. 

A  general  form  of  the  equivalent  circuit  representing  the  cavity 
and  its  coupling  network  is  shown  in  Figure  4.  The  cavity  system  con 
sists  of  a  cavity  and  part  of  the  input  transmission  line  to  which  the 
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Material  -.copper 
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cavity  is  coupled  by  a  loop.  To  be  complete,  the  equivalent  circuit 
should  contain  enough  elements  to  account  for  all  phenomena  asso¬ 
ciated  with  the  cavity.  Such  a  circuit  is  difficult  to  deal  with,  so  we 
first  construct  a  reasonably  complete  equivalent  circuit  and  then  at¬ 
tempt  to  simplify  it  by  making  physically  reasonable  assumptions. 

The  terminals  of  the  coupling  network  1-1  and  2-2  are  chosen  so 
that  all  sources  of  loss  are  included  between  them.  Also,  1-1  will 
be  located  at  the  detuned-short  position  (DSP)  in  order  to  simplify 
the  treatment  of  data.  Of  the  several  possible  forms  of  equivalent 
circuits,  the  one  shown  in  Figure  5  is  used,  since  it  is  most  readily 
analyzed. 

In  this  network,  the  resistances  R  and  R  represent  the 

s  p 

losses;  the  reactances  X  and  X  represent  the  reactive  elements, 

p  s 

such  as  the  inductance  of  the  coupling  loop.  A  simple  analysis  near 
resonance  is  possible  only  if  the  energy  storage  in  the  loop  circuit  is 
negligible  compared  with  that  in  the  cavity  at  resonance.  It  is  then 
possible  to  neglect  the  reactances  X^  and  X^  and  use  the  equivalent 
circuit  in  Figure  6,  which  has  a  new  resonant  impedance  of 
after  transformation  through  the  transformer. 

A  plot  of  the  variation  of  the  impedance  at  terminals  1  -1 
(plane  P^)  with  frequency  is  a  circle  on  the  Smith  Chart,  since  it 
represents  the  conformal  transformation  of  the  circular  impedance 
locus  of  the  network  at  plane  b-b  through  the  coupling  network.  The 
locus  is  shown  in  Figure  7.  If  the  series  resistance  in  Figure  6  were 
zero,  the  impedance  locus  would  pass  through  the  origin  (r  -  0  , 
x  =  0);  thus  the  choice  of  corresponds  to  the  location  of  the 
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detuned-short  position  in  the  case  of  no  loss. 

Referring  to  Figure  7,  let  the  intercepts  of  the  impedance 
locus  with  the  resistive  axis  be  called  a  and  p  ,  which  are  defined 
as 


a 


(9) 


where 


P  =  a  +  p 


P  = 


R  p, 

P  ^2 

R  +  p,Z 
p  rZ  c 


+  1 


(10) 


(11) 


y 


(1Z) 


The  intercept  described  by  a  results  when  the  cavity  is  tuned  far  off 


resonance  and  that  described  by  P  when  it  is  tuned  to  cj  ,  where  the 

o 

resonant  frequency  o>q  is  taken  as  the  frequency  when  the  impedance 
across  R^  is  a  maximum  and  purely  resistive. 

The  input  impedance  L ^  at  plane  P1  may  be  written  as 


R  + 


R 

p  bb 

R  +  Z,, 
p  bb 


(13) 


and 


Pz 

1  +  T~Z  0o  6 


(14) 
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v/h<-  re 


Z  -  chara.  !•  rs!"  muicdaiu  of  t  rafismiismn  I  i  i  •  • 
<> 

Q  -  an !o.i(i*-<i  O  <>t  vavitv  and 

n 

6  =((*.-(•»  ) ! w  -  tunino  parameter 

o 

Combining  Kqu-il:uns  (I  J)  ami  (\4)  and  rea  r  ranging  gi\e*~ 


tn  know  the  >oitp!:n>,  network  pa  ratnete  r  s  R  a  ml  K  to  in-  able  to 

mtiTtiM  tin:  nnpedar.ee  data  The  values  ol  K  and  V  <  ouid  be 
•  s  |> 

dele  rininvd  experimentally  if  tit*--  coupling  network  could  be  piiystt  1 :  i  y 
sepa rated  front  lit*-  vavity.  This  >s  impractical  m  i  to  pus  si  bio  iri  most 
cases.  Since  iti  most  devices  coupling  losses  are  **tnaii  and  arc  (lac 


mainly  So  either  ^-r:cs  *ir  shunt  losses  alone,  it  is  <  onvenient  to  ,  uii- 

4 

sider  three  separate  cases 
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A.  CASK  1.  SHUNT  LOSSES 


In  this  case  a  =  0  and  Equation  (15)  becomes 


1  +  j  2  Q  5 
J  o 


where  O  may  be  found  from  the  measured  impedance  locus;  and  the 
true  value  of  Q  is  obtained  from  Equation  (16)  if  y  is  known. 


It .  CASE  2.  SERIES  LOSSES 

In  this  case,  R  =  oo  ,  and  Equation  (15)  becomes 


II-  j  Z  U  b 
1  o 


The  value  of  O  may  be  obtained  from  the  inrastiml  imped,' mce  locus 
in  llic  following  way.  Let  2Q  f>  I  and  substitute  in  Equation  (18) 

to  get 


1  (1  I  a  J_  j  (|i  -  a) 


Using  values  of  <i  and  [1  taken  from  a  plot  ol  the  data,  one  can  com¬ 
pute  (he  real  and  imaginary  pails  ol  Equation  (IV)  and  locate  them  on 
the  locus  determining  the  points  for  which  2Q  6  =  _t  I  II  I  he  fre¬ 
quencies  at  which  they  incur  are  f  and  f  ,  then 


“l  -  f2 
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If  the  two  points  on  the  locus  do  not  coincide  with  data  points  for 
which  the  frequencies  are  known,  the  desired  frequencies  may  be 
found  by  constructing  a  linear  frequency  scale  as  shown  in  Figure  9. 


C.  CASE  3.  SERIES  AND  SHUNT  LOSSES 

When  both  series  .tnd  shunt  losses  arc  ore  sent  and  each  is 

t 

small,  the  Qq  value  may  be  found  from  Equation  (20),  which  takes 
into  account  the  series  losses  alone,  and  then  multiplying  this  value 
by  the  factor  given  in  Equation  (16)  to  account  for  shunt  losses. 
When  the  behavior  of  resonant  cavities  is  considered, 

(the  external  Q)  and  (the  loaded  Q)  are  also  of  interest. 


1  C 

These  along  with  Q  and  Q  are  defined  asJ 


energy  stored  in  cavity 
energy  dissipated  in  cavity 


'  _  energy  stored 

r  Q  fir  lincepC  4-  c 


cavity  losses  +  shunt  losses 


q  -  w  energy  stored  in  cavity 

cavity  losses  t  coupling  loss 
+  energy  dissipated  in  load 


Q  =  oj 


energy  stored  in  cavity 
energy  dissipated  in  load 


When  coupling  losses  are  not  present,  or  if  Z  ,  R  ,  and  R  are 

o  s  p 

considered  to  be  the  load,  Q  ,  Q.  ,  and  Q  are  related  by 

X  Li  O  1 


1  _  1  ,  1 

G7  "  W  +  O' 

-Li  OX 
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When  7.  and  R  are  the  load,  then 
o  s 


i  +  J_ 

^  Qx 


If  and/or  Qq  have  been  experimentally  obtained,  and  a  ,  p  ,  and  y 

are  known,  then  the  expressions  for  Q,  and  Q  are  found  as  follows. 

1  -L  x 

Let  V  be  the  voltage  across  the  capacitance  shown  in  Figure  6,  then 


from  Fc|ualion  (21), 


7  CV2 

«o  "  j  j— V2T  j  =  w  C  P?.Zo 


2  LP2Zo1 


Similarly,  from  (-equation  (lf>) 


O  U  — f-T  =  oiCp,7,  - 

o  o  y  I  1  i  o  y  +  1 


o'  LJ.  5  ,  Q _ 1 _ 11  ii 

c>  I  I  p  o  \  I  l  1  l  (3 


.nit]  1  i  t  mi  i,.(jii  at  i<  >i  i  (  /.()) , 


r  '  1  I  a 

Q  ■:  U  n - - 

x  o  p  -  a 


Thus  it  a  ,  p  ,  ami  y  n  re  all  known,  it  is  possible  to  ch*t»*  rm int*  tin* 
various  Q  values.  II  y  is  not  known,  Iho  Q  values  may  only  be 
lound  in  Irrms  ol  Q  ,  in  whirh  tin*  cavity  lo.ssos  and  sbunl  losses  in 


tin*  f  oiipling  network  art*  not  srpa  rail'd . 
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For  the  particular  case  of  interest,  it  is  convenient  and  rea¬ 
sonable  to  assume  that  the  shunt  losses  are  very  small  and  to  include 

I 

them  with  the  cavity  losses.  Then  from  Equation  (16),  Qq  and  Qq 
become  identical,  and  only  the  effect  of  series  loss  is  explicitly  con¬ 
sidered.  A  method  for  treating  the  data  is  given  in  Case  Z,  and  an 
alternate  method  follows.  Tne  data  required  are  a  plot  of  VSWR 
versus  frequency  near  resonance,  and  the  values  of  a  and  (3  .  The 
VSWR  at  the  half-power  frequencies,  where  2Qo&  =  +.  1  ,  is  given 

by 


I  1  +  j  +  (P  +  ia)  I  +  I  1  +  j  -  (P  +  jo-)  I 

I I  +  j  +  t(3  +T°Ti  1  1 1  +  j  -  (P  +  jHTT 


which  may  be  simplified  to 


(31) 


Z  \  (pZ  +  a*'}  I  J  4  I-  (p4  +  u4)  -  Hup  -I 

'V;  "  2(p  I  a) 

Z  u 


The  f rnqucm: Les  and  may  bo  found  from  the  plot  of  VSWR  vi-  rsus 

frequency,  and  .liquation  (ZO)  can  be  used  to  find  C2 

o 


IV.  MEASUREMENT  OF  i/  R  ,  /Q 

1  SO  o 

A  prediction  of  the  performance*  of  cavitiej s  its  resonant  systems 

is  made  possible  by  tne  measurement  of  R^  ,  the  shunt  resistance, 

and  Q  of  the  resonator,  or  by  the  measurement  of  R  ,  /Q  in 
o  HI1  o 
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the  case  where  an  electron  beam  is  present.  These  three  parameters 

are  considered  because  cavity  resonators  are  subject  to  three  types 

of  loading:  resistive  losses,  external  loads,  and  beam  loading.  For 

maximum  output,  when  the  cavity  is  used  to  couple  energy  from  an 

electron  beam  to  a  load,  it  is  customary  to  design  the  resonator  so 

2 

that  the  effective  shunt  resistance  p  is  equal  to  the  beam- 

loading  impedance;  which  may  be  approximated  by  the  d-c  beam 

impedance.  This  specifies  R  ,  and,  since  R  ,  /Q  is  inde- 

sh  sh  o 

». 

pendent  of  loading,  determines  Lhe  degree  of  coupling  necessary  for 

proper  matching  Appropriate  matching  will  be  necessary  to  obtain 

reasonable  output  power  from  the  beam -plasma  system,  therefore 

knowledge  of  p"  R^/Ci  for  the  cavities  used  is  useful. 

The  effective  gain -bandwidth  parameter  p^  R^  /Q  for  a 

cylindrical  electron  beam  of  radius  b  traversing  an  interaction  gap 

of  length  cl  in  a  resonator  having  cylindrical  symmetry  can  be  meas- 

5 

u  red  by  a  perturbation  technique.  The  resonant  frequency  shill 
caused  by  a  small  metallic  sphere  located  successively  aL  closely 
spaced  points  on  the  axis  of  the  interaction  gap  is  measured,  and  the 
resulting  distribution  of  frequency  shifts  is  used  to  determine  both 
the  configuration  of  the  elcitric  field  on  the  axis  of  the  resonator  and 
the  value  of  p*"  R  ,  /Q  that  corresponds  to  the  electron  beam  to  be 
used  in  the  resonator. 

The  equation  relating  p^  R  /Q  to  the  frequency  shifts  caused 

by  a  metallic  bead  of  radius  r  is 
'  o 
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where  u  is  the  stored  energy,  and  W  is  the  energy  dissipated  per 
second  :n  the  cavity  The  effective  shunt  resistance  is  defined  in 
terms  of  the  square  of  the  effective  gap  voltage  as 


R 


eil 


i£ 


ffective  voltage)2 
power  loss 


(35) 


> 


{^b  ^d  Em  d) 

— m - 


(36) 


where  and  ^  are  as  defined  by  Branch.  Combining  Equations 

(34)  and  (io)  gt\es 


!.  „  .  '“b  “d  S,,'11' 

P  K  .  /Q  - 
s  h  o 


Tznr 


(37) 


i  . 


From  Staler ‘s  perturbation  formula  for  a  metallic  sphere, 


3  f(n  H2  -  c  E2)  dV 
-  t  i  r  *• 

o  J 

7 -  =  - TD - 

o 


(38) 


and  for  the  re-entrant  cavities  of  interest,  H  =  0  on  the  axis.  Since 
the  head  is  small.  Equation  (38)  then  becomes 


Af  3  e  E2  AV 

T  = - TV - 

o 


(39) 


Solving  Equation  (39)  for  U  and  substituting  in  Equation  (37)  gives 
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(40) 


where  g(z)  is  the  normalized  variation  of  the  field  in  the  gap  referred 
to  the  field  on  the  axis.  Thus  for  an  arbitrary  field  variation 

E 

gU)  =  jp-  ■  (45) 

m 


From  Slater's  perturbation  formula,  the  change  in  frequency  at  any 
point  on  the  z-axis  is  proportional  to  the  square  of  the  electric  field; 
thus 


E  =  K*/Af 

z 


E 

m 


K  \TaT 


max 


g(z)  = 


nTa? 


max 


(46) 


Substituting  gives  Equation  (33): 


2  Rsh 

^  TJ- 

o 


=  120  x 10 


8  o 


i2  (*M 


i2(V>) 


f 

-0O 


_  jpz 

*/Sf  e  dz 


ST 


(47) 


The  integral  can  be  evaluated  graphically  or  by  Simpson’s  rule, 

A  diagram  of  a  test  setup  suitable  for  performing  the  required 
measurements  is  shown  in  Figure  10.  The  resonance  curve  is  dis¬ 
played  on  the  oscilloscope;  then  the  marker  signal,  an  accurately  known 
harmonic,  is  centered,  and  the  frequency  is  read  on  the  electronic 
counter.  The  bead  is  then  lowered,  the  marker  signal  is  readjusted 
to  center,  and  the  corresponding  frequency  shift  is  read  from  the 
counter.  The  process  is  repeated  in  convenient  steps. 
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est  S»*tup  for  Perturbation 


The  equations  and  techniques  developed  in  this  section  may  be 

used  to  find  R  ,  /Q  for  the  cavity  alone,  as  well  as  p2  R  .  /CJ  when 
sh  o  8,11  o 

an  electron  beam  must  be  considered. 


V.  COLD  TEST  RESULTS 


The  first  measurements  were  made  on  a  test  cavity  built  for 

cold  testing  and  the  results  are  compared  with  theoretical  values. 

Q  was  measured  using  a  coaxial  slotted  line  and  a  standing-wave 
o 

indicator,  the  initial  value,  obtained  by  following  the  procedure  dis¬ 
cussed  in  Section  III.  being  Qq  =  4920.  Lower  values  were  obtained 
from  later  measurements  because  of  the  deterioration  of  the  copper 

surfaces.  R  /Q  was  found  by  using  the  perturbation  technique 
sn  o 

described  in  Section  IV.  Dividing  both  sides  of  Equation  (41)  by  p 
gives 


120  x  10 


8  d^  ^max 

~Z  AV"“ 


(48) 


A  metallic  bead  was  pulled  through  the  gap  and  the  maximum  frequency 
shift  recorded.  The  following  values. 


d 

-  0.125  in.  -  0.318cm 

AV  =  4irr  /  3 
o 

f 

o 

=  3280.97  Mc/s 

r  =  0.031  in.  =  0.0787  cm 
o 

Af 

=  94  kc/s 

max 
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inserted  in  Equation  (48)  gave  a  value  of  51  for  R  j  /Q  •  For  com¬ 
parison,  theoretical  values  follow,  together  with  the  values  for 

g 

O  and  R  ,  /Q  obtained  from  Ginz.ton  and  Nalos: 
o  sh  o 


f  (Mc/s) 
o 

o! 

a  | 

R  ,  /Q 

sh  o 

Theoretical  Value: 

Experimental  Value: 

3260 

3280. 97 

5250 

49  20 

45 

51 

The  close  agreement  between  the  theoretical  and  experimental  values  is 
an  indication  that  the  design  procedures  were  sound. 

Figure  11  shows  the  normal i/.ed  axial  electric  field  in  the  gap 
as  determined  by  the  perturbation  method.  The  curve  is  not  symmetri¬ 
cal  because  the  gap  configuration  of  the  singly  re-entrant  cavity  is 
asymmetrical. 


VI.  FREQUENCY  SHIFT  OF  CAVITY  MODES 
IN  PRESENCE  OF  A  PLASMA 

In  this  section,  an  expression  for  the  frequency  shill  of  a 
nondogenorate  mode  of  a  ravily  perturbed  by  a  plasma  is  derived. 

The  properties  of  the  plasma  are  represented  by  a  tensor  dielectric 
constant  or  by  a  tensor  conductivity  derived  from  the  following  tensor 

mobility:  ^ 
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Cavity  orientation 


V  +  r 


-J(1  -  r) 


0 


!  j(t  -  r) 


()  V  _ 

v  +  j(u)  -  to. 


/  3  Of  . 

±i  /  v  ^  ' 

^  j  vc  ^  j(w  - 


V  -i  a  t  . 
in  j  V  JH-  ^ 
.i  I  V  i  jui 


,iih1  >  11  /in  in  tin-  electron  ey<  lot  rmi  frequency;  B  is  the* 

(I  i  jn<i|>in'l  it-  field  parallel  to  tlie  -/,-ax.iK,  and  f  is  Llie  normalized, 
.spin-  rii  al  ly  sy  mi  net  ric  pa  it  ol  the  electron  velocity -distribution 
1 1 1  i  >  ■  t  ion.  Winn  f  in  Maxwellian  and  v  is  independent  of  velocity, 

V.. 

Mu  i  i  >1  fj; » ••  hi  fin*  .  i  ,  and  |>  tan  In*  evaluated  to  yield 


vi:  h  j(w  1  ^ 


v  -i  j[to  - 


V  I  joj 


The  mobility  and  conductivity  are  related  by 


ne  |i  = 


0r 

0 


1  +  r 


jtt  -  r) 


re 

°ee 

o 


-3(1  -  r) 


l  + 


2p 


(61) 


where  n  is  the  electron  density.  A  tensor  dielectric  constant  is 
defined  by 


K  -  1  i  a/jioe. 


K 


K 


Or 


K 


i-O 


K 


(10 


(62) 


where  I  is  the  unit  tensor, 


K  =  K 
CIO  r 


'•>'  0  -  jv' 


(1  -  jv) 
- ^ 


( (>  hi ) 


Kh!  =  -  K0r 


2,2  , 

Ki)  f  hi  -  >i).  f  (0 

p _ b _ 

J  (I  -  ,i/j  -  «■?/</ 


«•  ''!>) 
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1 


1 


(63c) 


2  2 

and  w  =  ne  /me  >  and  v  =  v  /w  .  Because  the  motion  of  ions  has 
p  o  c 

been  neglected,  the  validity  of  this  treatment  is  limited  to  plasmas  of 

electron  temperatures  low  enough  that  Landau  damping  and  other 

effects  characteristic  of  a  hot  plasma  are  negligible.^ 

We  now  consider  two  solutions  to  Maxwell's  equation  for  a  cavity: 

E  ,  H  in  the  absence  of  the  plasma,  and  E,  ,  H,  when  the  plasma 
o  o  11  1 

is  present.  The  first  set  satisfies 


V  x  K  =  -  jio  u.  fl 

O  0*0  0 


V  x  It  -  jw  e  (K  •  E  )  , 

o  o  o  =  o  o 


where  K()  is  the  unit  tensor  if  the  cavity  is  empty.  When  the  plasma 


is  present,  E,  ,  II.  satisfy 


V  x  Ej  =  -  jw  |io  Uj 

v  X  Tlj  -  (K  ■  Ej) 

Now  multiply  Equation  (65)  by  E ^  and  Equat  ion  {(>1)  by  E  ,  and 
integrate  over  the  cavity  volume  to  obtain 


V  x  If  dv  =  j 


j  I  w  c 

f  o  o 


E,  •  (K  .  E  )  dv 
1  =o  o 


E  ■  v  x  H  dv  e  j  i  oj£  E  (K  ■  E.)  dv 
o  I  J  o  o  =  -1 


Thin  subtrai  t  Equation  (G'l)  from  Equation  (6H),  which  give 
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J  [JE,  '  V  x  H  -  E  •  V  X  H.J  dv  =  j  ( €  E, 

J  1  O  O  1— ^  I  o  ^  o  1 


(K  •  E  ) 
o  o' 


-  u>  Eq  •  (K  •  ExOdv  .  (70) 


The  left  side  of  Equation  (70)  becomes 


[EJ  •  V  x  Hq  -  Eq  •  V  x  dv  =  J  V  C  Ho  x  ^  -  Hj  x  EqJ  dv 


+  |  CH  •  V  x  E.  -  H.  ■  VxE  ]dv 
j  o  1  1  o 


and.  by  the  divergence  theorem,  the  first  term  on  the  right  of  Equation  (71) 
can  be  rewritten  as  the  surface  integral  over  the  cavity  walls  of 

n  ■  [H  x  E.  -  B,  x  1  ]  , 

^  o  1  1  o  J 

which  is  zero  if  the  walls  are  perfect  conductors.  The  second  integral 

in  Equation  71  can  be  simplified  by  using  Equation  (04)  and  Equation  (66) 

to  define  and  vx  • 

o  1 


[  CH  •  VxE.  -  HT  VxE  ]dv  =  ju  (u  -  u)  f  H  ,*  H 
1  1  °  r°  o  J  1  o 


dv  .  (72) 


Equation  (70)  can  now  be  rewritten  as 


>.  (w  -u)H,  ■  H  dv  =  £  E,  ■  (K  •  E  )  -  u  E  •  (K  •  E,)l  dv 

O  O  'J  1  O  m  j  *—  O  1  =0  o  o  =  1' - 


which  can  be  solved  for  the  relative  frequency  shift  6  =  (oj  w  ) / to  as 


e  [E  ■  (K  '  E.)  -  E.  •  (K  ■  E  )]  dv 
o  —  1  1  — o  o 


(74) 


e  I  E.  ■  (K  '  E  )  dv  -  u.  H.  ■  H 
cO  i  l  =o  o  ro  /  1  o 


dv 


Thin  can  also  he  written  as 


6 


■  i_(tr  -  <r  ) 
=  -o 


0) 


dv 


Ej  J  dv 

-  u  H  •  H,  dv 
1  o  J  o  1 


(75) 


For  the  case  where  the  unperturbed  cavity  is  in  a  vacuum,  Equation 
(75)  can  he  further  simplified  to  give 


6 


j 

iu)c 

O 


I  E  •  n  •  E  dv 

J  o  =  1 


tW 


(76) 


This  equation  is  the  basis  of  the  plasma  diagnosis,  since  it  relates  the 
complex  frequency  shifts  to  the  cavity  fields  and  plasma  parameters. 


VII.  FI, ASM  A  DIAGNOSIS 

In  the  microwave -cavity  method  of  measuring  properties  of  a 
plasma  (in  use  for  over  fifteen  ycars^  ’  ^),  the  plasma  to  be  studied  is 


partially  or  completely  contained  in  a  microwave  cavity,  and  the  resulting 


shift  in  the  resonant  frequency  and  the  change  in  cavity  are  measured. 

These  two  measurements  are  then  related  to  the  electron  density  and 
electron  collision  frequency  by  means  of  the  expressions  derived  in  the 
preceding  section. 

The  structure  of  the  singly  re-entrant  cavity  makes  it  possible 
to  assume  that  only  z-directed  electric  fields  interact  strongly  with  the 
plasma.  This  means  that  Equation  (7b)  may  be  simplified  to 
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Ur  E  '  E.  dv 
I  zz  o  1 

J  V  E1  dv 


wh<*  re  is  given  by 


tv 
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iu;”'  1  __  ne  y  -  j 

~  Tj- T'jI3"  ~  ~"~2 


(77) 


(78) 


Sim  e  o  is  a  complex  quantity,  f>  is  also  complex,  with  its  real 

part  equal  to  the  relative  shift  in  resonant  frequency,  and  with  twice 

the  itnajpnaty  pin  I.  equal  to  i  .  The  following  equations  express 

<  > 

these  r  el  at  ionshipsl 
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1  X 


P.J  dv 
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2  2  2  2 

where  x  =  ne  /me  <o  =  /u  and  V  =  v  / cj  . 

o  p  c 

The  only  problem  remaining  is  that  of  performing  the  required 

1  3 

integrations  in  Equations  (79)  and  (80).  Buchsbaum  and  Brown  have 
solved  this  problem  by  assuming  that  the  plasma  causes  only  a  small 
perturbation  of  the  cavity  fields  and  the  electron  density  is  constant 
within  the  volume  of  integration.  These  assumptions  allow  one  to  re¬ 
place  E,  with  E  to  obtain 
r  1  o 


Af 

f 


(81) 


E2  dv 
o 


Itf 


2  2 

When  the  electron  density  is  low,  that  is,  when  x  =  cWu>  «  1  ,  the 

approximation  is  good,  since  the  plasma  does  not  appreciably  disturb 

the  electric  field.  As  the  density  is  increased  however,  Ej  becomes 

different  from  E^  for  three  reasons:  (1)  a-c  space  charge,  commonly 

2  2 

called  "plasma  resonance,"  appears  as  u>  / u>  approaches  unity;  (2) 

2  2 

as  to^/u)  ‘  becomes  greater  than  uniLy,  the  plasma  begins  to  shield  its 
interior  from  the  field  outside;  and  (3)  when  both  w2/u>2  and  the  pres¬ 
sure  arc  high  enough  that  the  Q  value  changes  greatly,  the  overlapping 
of  higher  modes  may  also  cause  to  be  different  from  E  by  adding 

to  E  some  of  the  fields  of  the  higher  modes.  This  means  that  the 
o 

cavity  method  will  work  well  at  low  pressures  and  small  percentages 
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of  ionization,  but  will  give  increasingly  poorer  results  as  pressure  and 
ionization  increase.  When  the  perturbation  theory  is  valid,  the  electron 
density  and  the  resonant  frequency  shift  have  an  approximately  linear 
relationship.  The  point  at  which  the  measured  data  begins  to  depart 
from  this  linear  relationship  indicates  the  beginning  of  the  region  where 
Equations  (81)  and  (82)  are  no  longer  valid.  Experimentally,  this  point 
obtained  by  making  a  plot  of  resonant  frequency  shift  versus  discharge 
cu r  rent . 

In  the  region  where  the  perturbation  analysis  is  not  valid,  one 
must  return  to  Equations  (79)  and  (80)  and  employ  the  actual  cavity 

fields  E  and  E,  if  an  exact  analysis  is  to  be  made.  This  is  a  ques- 

o  1 

tionable  procedure  because  the  approximations  that  have  to  be  made  in 
order  to  obtain  a  reasonable  solution  fur  the  fields  are  ■  lif f icu It  to  work 
with.  For  convenience  in  experimental  work,  it  may  be  preferable  to 
extrapolate  the  data  taken  in  the  region  where  the  perturbation  relations 
are  valid  and  then  use  the  point  where  u)  =  for  calibration.  The 
validity  of  this  method  can  only  be  determined  by  experiment. 

Since  the  electron  density  is  not  a  constant  in  the  actual  system, 
hut  varies  radially,  the  values  for  n  found  by  experiment  will  be  aver¬ 
age  values.  These  are  not  a  disadvantage  in  this  study  since  we  are 
concerned  with  the  bulk  plasma  properties. 
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VIII.  EXPERIMENTAL,  RESULTS 


Preliminary  data  has  been  obtained  on  the  hydrogen  plasma 
discharge  described  in  this  report.  Figures  12  and  13  show,  respec¬ 
tively,  the  d-c  current- voltage  characteristics  and  the  detuning  of 
the  microwave  cavity  as  a  function  of  discharge  current. 

The  linear  relationship  between  the  cavity  detuning  and  the  dis¬ 
charge  current  indicates  that  the  approximations  in  the  perturbation 
analysis  are  justified  for  the  range  of  discharge  currents  studied. 
Further  measurements  are  being  made  to  find  the  absolute  electron 
density. 

II.  was  possible  to  reset  the  plasma  current  so  that  the  cavity 
detuning  was  within  one  or  two  megacycles  of  the  earlier  data.  These 
results  indicate  that  cavity  detuning  is  a  practical  way  of  adjusting 
the  plasma  density  to  the  desired  point  lor  beam  plasma  amplifications. 
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ABSTRACT 


The  results  of  measurements  of  the  properties  of  trans¬ 
verse  and  radial  ion-induced  electron-beam  oscillations  are 
described  in  detail.  From  the  data  on  the  transverse  oscillations 
and  from  data  published  by  Mihran  and  Agdur,  it  is  concluded  that 
the  transverse  oscillations  increase  in  amplitude  with  time 
because  of  interaction  ajter  the.  electrons  enter  the  magnetic 
focusing  field.  Because  the  ion  core  disturbs  the  potential  well 
caused  by  the  electrons,  the  electron  beam  follows  a  helical  path 
from  the  gun  to  the  collector.  Ions  or  electrons  traveliiig  toward 
the  gun  alter  the  helical  path  and  cause  the  oscillations  to  grow. 

The  radial  oscillations  were  found  to  occur  at  a  frequency 
th;iL  was  1.4  to  2  times  the  transverse  frequency.  They  increase 
in  amplitude  with  time  because,  accompanying  the  oscillation  of 
the  ion  core,  there  is  a  density  modulation  of  the  ions  in  the 
core  that  affects  the  focusing  of  the  electron  beam  in  the  conver¬ 
gent  Pierce  gun  in  a  periodic  manner.  This  causes  the  diameter 
of  the  electron  beam  to  oscillate,  which  in  turn  causes  the  radial 
oscillations  of  the  ions  to  grow. 

Suggestions  are  made  for  the  control  or  elimination  of 
radial  and  transverse  oscillations. 
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I.  INTRODUCTION 


Although  a  large  amount  of  analytical  and  experimental  work 

has  been  carried  out  during  the  past  35  years  in  an  attempt  to 

explain  ion-induced  electron-beam  oscillations,  the  questions  of  how 

and  why  the  oscillations  build  up  remain  unanswered.  The  purpose 

of  this  paper  is  to  present  the  results  of  detailed  experimental 

analyses  of  the  behavior  of  the  electron  beam  during  oscillations. 
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These  results  verify  and  complement  the  work  of  Mihran  and 
Agdur,1^  which  has  provided  the  greatest  amount  of  experimental 
data  to  date  on  ion-induced  electron-beam  oscillations. 

In  addition,  an  extension  of  the  analytical  work  of  Peterson 
1 1 

and  Puthoff  is  given,  which  helps  in  explaining  the  build-up  of 
oscillations.  Peterson  and  Puthoff  treat  the  transverse  and 
radial  oscillations  of  a  cylindrical  ion  cure  within  a  cylindrical 
electron  beam  that  is  assumed  to  be  stationary  (the  time  spent  by 
electrons  in  tiie  plasma  is  assumed  to  be  very  short).  The  frequency 
of  transverse  oscillations  of  the  ion  core  is  found  to  be 


where  <u  is  the  plasma  ion  frequency  first  found  by  Langmuir  and 

1 

Tonks,  c  is  the  radius  of  the  ion  core,  and  d  is  the  radius  of  the 
anode  within  which  the  oscillations  are  occurring.  The  frequency  of 


1 


radial  <:  sri J  latii .u  of  the  ion  core  is  found  to  be  simply  ,j,,  when  the 
amplitude  of  the  oscillations  is  small. 

From  the  analytical  and  experimental  results,  the  mechanism 
for  the  bnild-up  of  ion  oscillations  is  deduced,  and  precautions  that 
maybe  taken  to  prevent  the  build-up  are  pointed  out. 


f.  - 


II.  DESCRIPTION  OF  EQUIPMENT 

The  beam  analyzer  used  in  obtaining  the  data  presented  in 
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this  study  is  the  one  described  by  Gilmour  and  Hailock  .  The 
electron  hcatn  on  which  measurements  were  made  was  generated  by 
the  shielded  pierce  gun  shown  at  the  right-hand  side  of  Figure  1. 

This  gun  was  normally  pulse  operated  at  S000  volts  and  had  a  micro- 
perveance  ol  1.  PS.  At  pressures  in  the  range  of  10  8  Torr,  the  area 
compression  of  the  beam  was  SO:!,  and  the  magnetic  field  required 
for  lirillouin  tlow  (less  Ilian  i.  per  cent  scalloping)  was  2.85  gauss. 

The  parts  of  the  analyzer  essential  to  the  current  measure¬ 
ments  are  shown  in  Figure'  1.  A  ball,  valve  is  shown  positioned 
between  I  lie  gun  and  the  drift  region,  which  contains  a  doubly  re¬ 
entrant  cavity  and  a  beam-ana  I  yzing  mechanism.  This  valve  is 
considered  to  lie  one  of  the  most  important  parts  of  the  analyzer 
because  when  it  is  cloned,  the  drill  tube  can  he  opened  to  atmospheric 
pressure  while  the  rul  Inule  remains  evacuated  at  a  pressure  of  about 
10  Torr,  Aftri  l  In-  drift  luhe  v/.i  s  rviinial.ctl,  the*  valve  was  opened, 
,mk1  the  Ihmmi  f  mm  I  lie  gun  passed  through  the  aperture  in  the  ball  to 
the  be.un-niochil.i  t  i  ng  and  bea  tn -ana  1  yzing  region. 

Tlie  m-H  ne-fi  n«».l  yzing  mechanism  shown  at  the  left  of  Figure  1 
is  shown  in  c  ross  seel  ion  in  Figure*  £.  Most  of  the  electron  beam  was 
collected  by  the  beam -c  ol  iecting  plate,  which  was  carbonized  to  reduce 
the  sec  wither y  emission.  A  small  portion  however,  passed  through  the 
.010 -inch  aperture  in  the  renter  ,»f  the  j>’  ite  to  a  Faraday  cage.  The 
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FIGURE  1.  Simplified  Sketch  oi  Beam  Analy 


cage,  in  turn,  was  attached  to  the  end  of  the  center  conductor  of  a 
nonresonant  coaxial  line,  which  was  used  to  measure  currents  at 


X-Y  recorder,  so  that  plotting  of  current  ampiitudeas  a  function  of 
position  in  the  beam  was  automatic. 

To  obtain  plots  of  the  current  density  versus  position  at 
various  times  after  the  beginning  of  Die  cathode  pulse  (necessary 
because  ions  caused  the  beam  shape  and  position  to  change  with 
time),  a  Hewlett  Packard  oscilloscope  with  a  display  scanner  was 
used.  The  output  from  the  display  scanner  was  applied  to  the  Y 
input  of  the  recorder. 

Pressures  up  to  approximately  10  ^  Torr  were  obtained 
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by  first  pumping  the  analyzer  to  the  low  range  of  10  Torr  with 
the  IS  l/s  ion  pump  and  then  bac  kfilling  to  the  desired  pressure 
with  hydrogen  through  a  palladium  leak.  A.t  pressures  in  the  range 
of  10  and  low  10  ^  Torr  the  "IS  1/s  ion  pump  was  used  to  help 
control  the  pressure;  at  higher  pressures,  the  1  l/s  ion  pump 
attached  to  the  election  gun  was  user!. 


III.  EXPERIMENTAL  OBSERVATIONS 


A.  GENERAL 

At  low  pressures  (  <  10  ^  Torr)  the  electron  beam  used  for 
this  study  was  very  nearly  an  ideal  Brillouin  beam  (at  a  magnetic 
field  within  Z  per  cent  of  the  Brillouin  value,  less  than  Z  per  cent 
scalloping  occurred).  No  ion-induced  oscillations  and  no  phase  shift 
of  a  r-f  signal  applied  to  the  beam  were  observed  during  cathode 
pulses  with  durations  up  to  30  ps.  This,  of  course,  indicates  that  at 

_7 

pressures  below  10  mm  Hg,  the  number  of  ions  produced  was 
negligible.  As  hydrogen  was  allowed  to  leak  into  the  analyzer,  the 

effects  of  ions  were  first  observed  at  a  pressure  of  approximately 
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3  x  10  'l'orr,  when  it  was  noted  that  the  beam  diameter  decreased 

IS 

slightly  as  a  luncLion  ol  time.  As  the  pressure  was  increased 
further,  the  beam  diameter  continued  to  decrease  until,  a  I.  a  pressure 
near  Z  x  10  Torr  a  few  cycles  of  small -amplitude  oscillations  with 
a  frequency  of  6  Mc/s  occurred  at  the  end  of  the  30-ps  cathode  pulses. 
These  osi  illations  were  of  the  t  ransverse  type  and  are  described  m 
detail  in  the  following  pages.  As  the  pressure  was  further  increased, 
the  transverse  oscillations  occurred  earlier  in  the  pulse  and  usually 
increased  in  amplitude.  These  oscillations  were  usually  stable  and 

were  amplitude  modulated  at  frequencies  up  to  approximately  Z  Me,  s. 

-4 

At  pressures  in  the  low  range  of  10  Torr  the  oscillations 
were  still  mostly  transverse,  with  a  few  cycles  of  radial  oscillat  ions 


Torr 


occurring  early  in  the  pulse.  At  a  pressure  of  5  x  10 
and  a  time  of  20  ps  after  the  beginning  of  the  pulse,  an  abrupt 
change  occurred  in  the  oscillations-they  changed  from  transverse 
to  radial  oscillations.  Usually  the  radial  oscillations  would 
remain  stable  for  a  few  cycles  and  then  become  very  unstable. 

At  the  same  time  that  the  radial  oscillations  became  unstable, 
strong  electron  cyclotron  oscillations  occurred.  The  cyclotron 
oscillations  are  not  described  in  this  report. 

B.  TRANSVERSE  OSCILLATIONS 

Figure  4  shows  records  of  the  electron  current  entering 

the  Faraday  cage  when  the  aperture  in  the  scanning  mechanism 

was  placed  at  seven  positions  on  a  diameter  of  the  beam.  Note 

the  phase  shift  of  180'  between  the  curves  at  the  top  and  those 

at  the  bottom  of  the  figure.  These  curves  and  the  curves  shown 

in  Figure  5  of  current  density  as  a  function  of  radial  position  for 

the  times  t ^  through  t^.  indicated  on  Figure  4  show  that  the  electron 

beam  was  moving  transversely  during  the  oscillations.  These  data 
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were  taken  at  a  pressure  of  4  x  10  Torr,  at  a  distance  of  15  inches 
from  the  electron  gun  with  the  beam  voltage  set  at  S  kv  and  the 
magnetic  field  at  320  gauss. 

Figure  6  shows  a  curve  relating  the  frequency  of  the 
transverse  oscillation  to  the  electron  accelerating  potential.  Note 
that  the  slope  of  the  curve  is  nearly  one-half,  which  is  the  value 
expected  for  ion-plasma  oscillations  when  the  ion  density  is 
pr  i portions  1  to  the  density  of  the  electrons  in  the  beam  from  the 
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electron  gun.  At  one  point  during  measurements,  the  hydrogen 
in  the  vacuum  system  was  replaced  by  nitrogen,  and  it  was  noted 
that  the  frequency  decreased  from  5.  6  Mc/s  to  1.5  Mc/s,  The 
ratio  of  these  two  frequencies  is  3.  15,  which  is  the  square  root 
of  the  ratio  of  the  mass  of  to  the  mass  of  and  which  is 
also  Lo  be  expected  for  ion-plasma  oscillations. 

Measurements  were  made  to  find  out  if  the  oscillation 
frequency  was  affected  by  the  strength  of  the  magnetic  field. 

The  unexpected  result  shown  in  Figure  7  is  related  to 
the  dependence  of  the  diameter  cf  the  electron  beam  on  the 
magnetic  field.  At  low  levels  of  magnetic  field,  the  beam  di¬ 
ameter  was  inversely  proportional  to  the  field,  whereas  at  high 
levels  (500  gauss  and  above)  the  beam  diameter  was  nearly 
constant  because  a  small  amount  of  magnetic  field  linked  the 
cathode^.  Since  the  beam  diameter  was  a  function  of  the  magnetic 
field,  the  shape  of  the  potential  well  and  the  ion-oscillation 
frequency  were  a  function  uf  the  magnetic:  field.  When  a  small 
coil  was  placed  around  the  electron  gun  and  was  used  to  intro¬ 
duce  magnetic  flux  into  the  cathode  region,  it  was  possible  to 
change:  the  shape  of  the  electron  beam  so  that  the  frequency  of 
oscillation  was  as  low  as  3  Mc/s  or  as  high  as  10  Mc/s,  while 
the  beam  voltage,  the  magnetic  focusing  field,  and  the:  pressure 
remained  constant.  It  was  also  possible  to  change  the  frequency 
of  oscillation  simply  by  changing  the  phase  of  the  cathode  pulse 
with  respect  to  the  a-c  filament  current  (normally  the  cathode 
was  pulsed  at  f>0  c.ps  with  the  pulses  being  phased  so  that  they 
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BEAM  VOLTAGE  (KV) 


FIGURE  £>.  Transveise  Oscillation 
Frequency  as  Function  of  Beam  Voltage. 


FIGURE  7.  Transverse  Oscillation 
Frequency  as  Function  of  Magnetic  Field. 
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occurred  at  zerois  •&:  the  =sn.-as'ai&<L» S  it i-iracut  current.)  Thus 
there  was  XiO  perl  urination!  <o< B  lihe  beam  by  the  magnetic  field 

OI  t he  fllci ucuX  . 

Measurements  mere  mnwde  «vi  the  frequency  of  oscilla¬ 
tion  a  a  a  function  01  jirc^izric  amid  time  at  a  beam  voltage  of 
5  kv  sit d  a  magnetic  t'leiliO  mi J(Q)  gauss.  The  frequency  at  the 
beginning  oi  transverse*  «r»«i«c u S D^i (t ii <:i>au»  was  found  to  be  constant 
at  8  Me/*  over  the  preii.»uie  range  of  S  x  10  —  S  x  10  Torr. 

•Ml.  <V 

A  marked  var.ail:  tn  an  d  re^iiine-mM'  \i  •  u  cur  red  as  a  function  of  time 
aft  er  the  be.j’sjvnw#  '•>»  ose  liUilatliii'  W*  (see  Future  ft). 


Fi-GCKE  «  Frequency  oi  T  raas  verse  Oscillations 
as  Funvtaoji  o:  T sane  «*tter  Beginning  of  Oscillations. 


Since  it  seemed  unlikely  that  the  beam  was  moving 
back  and  forth  in  a  straight  line  during  transverse  oscilla¬ 
tions,  complete  plots  of  the  current  density  as  a  function  of 
horizontal  and  vertical  position  were  made  at  60-degree 
intervals  during  a  few  cycles  of  the  transverse  oscillations. 
These  plots  were  made  near  the  beginning,  of  the  current 
pulse,  so  that  the  build-up  of  the  transverse  oscillation 
could  be  observed.  The  center  position  of  the  beam  was 
determined  as  a  function  of  time  from  the  plots  of  current 
density  (see  Figure  91-  ffColtce  that  the  beam  first  moved 
in  a  path  that  was  nearly  linear,  but  with  time,  the  oscilla¬ 
tion  path  became  elliptical.  (The  data  for  Figure  9  was 
taken  at  an  axial  position  4  in.  from  the  electron  gun.  ) 

With  the  data  from  Figures  4.  5.  and  9;  it  is  clear 
that  at  any  axial  position,  the  electron  beam  followed  an 
elliptical  path  during  transverse  oscillations,  ft  is  interesting 
to  observe  that  because  of  the  low  transit  time  of  electrons 
from  the  gun  to  the  collector  (approximately  .  10  x  10  ^s),  no 
single  electron  moved  in  the  elliptical  orbit  of  the  brain, 

because  the  time  required  to  traverse  one  orbit  was  about 
-6 

.15  x  10  s;  thus  the  assumption  of  a  stationary  beam  by 

11  .  _ 

Peterson  and  Puthoff  is  at  least  partially  justified.  To 

determine  just  how"  stationary"  the  beam  was  as  it  moved 
through  the  beam -analyzing  regions,  current  density  was  plotte 
as  a  function  of  radial  position  for  various,  axial  positions  from 


the  gun  at  times  corresponding  to  t^  and  (see  Figure  4). 

From  these  plots,  the  center  position  of  the  beam  in  a  hori¬ 
zontal  plane  as  a  function  of  axial  position'  was  determined  and 
is  shown  in  Figure  10.  If  the  cyclotron  oscillation  in  Figure  10 
is  ignored  it  is  apparent  that  instead  of  following  a  linear  path, 
the  beam  followed  a  path  that  was  a  portion  of  a  relatively  low- 
frequency  {compared  with  the  cyclotron  frequency)  sinusoid. 

Figure  11  shows  this  path  more  clearly  and  shows  that  it  was 
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actually  helical.  This  result  has  also  been  observed  by  Agdur. 

The  reason  for  the  path  being  helical  instead  of  linear  may  be 
understood  if  the  effect  on  the  electron  trajectories  of  an  electron- 
beam  well,  which  is  nearly  filled  with  ions,  is  examined.  Let 
us  assume  that  the  ion  core  is  displaced  to  one  side  and  the 
electron  beam  to  the  other  side  of  the  axis  of  the  solenoid,  as 
is  shown  in  Figure  12a.  Then  a  very  rough  approximation  to 
the  potential  distribution  in  the  electron  beam  as  a  function  of 
radial  position  from  the  axis  of  the  solenoid  is  the  parabolic 
distribution  shown  in  Figure  12b  ,  where  it  has  been  assumed 
that  neutralization  has  occurred  in  the  ion  core.  If  it  is  assumed 
that  no  magnetic  flux  links  the  cathode,  then  it  may  easily  be 
shown  that,  the  trajectory  of  an  electron  is  of  the  form, 

r  =  A  sin  (0  1  Ct)  (1) 

where  0  is  t  (co^  is  the  Larmor  frequency)  and  C  is  a  constant 
determined  by  the  density  of  the  beam.  Equation  { 1 )  shows  that 
the  electron's  path  is  a  circle  which  precesses  about  the  solenoid 


17  - 


axis.  If  the  constant  C  is  evaluated  for  a  potential  drop  of  10 
volts  between  the  center  of  the  solenoid  and  the  outer  edge  of 
the  beam,  the  precession  rate  shown  in  Figure  10  is  obtained. 
This  10-volt  drop  is  quite  conceivable,  when  it  is  remembered 
that  the  potential  drop  from  the  edge  to  the  center  of  the  beam 
considered,  with  no  ions  present,  was  70  volts.  From  the. 
preceding  considerations,  it  is  apparent  that  the  transverse 
oscillations  originated  before  the  electron  beam  entered  the 
region  in  which  beam  measurements  were  made. 

It  had  been  shown^  tha.t  under  some  circumstances 
secondary  electrons  from  the  electron  beam  collector  provide 
the  feed-back  loop  necessary  for  the  build-up  of  transverse 
oscillations.  For  the  measurements  made  in  this  study,  this 
occurred  at  low  pressure  (<  5  x  10  Torr).  As  the  collector 
was  moved  transversely  or  axially,  both  the  amplitude  of  the 
oscillations  and  the  time  at  which  they  started  were  greatly 
affected,  even  though  the  collector  was  coated  with  Aquadug 
to  reduce  the  emission  of  secondary  electrons  to  a  minimum. 

It  seems  probable  that  it  was  the  high-speed  rather  than  the  low- 

speed  secondary  electrons  that  were  responsible  for  the  feed- 
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back  action,  because  in  the  pressure  range  <  5  x  10  Torr  for 
microsecond  pulses,  the  beam  was  not  completely  neutralized, 
causing  a  potential  depression  as  high  as  a  few  hundred  volts 
in  the  drift  tube.  This,  of  course,  was  enough  to  stop  the  low- 
speed  secondary  electrons.  By  using  a  special  beam  collector, 
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Agdur  had  trapped  high-speed  s.s  well  as  low-spccd  secondary 
electrons,  but  although  he  succeeded  in  reducing  the  amplitude 
of  the  transverse  oscillations  by  a  factor  of  5,  he  was  not  able 
to  eliminate  them. 

Mihran  found  that  the  suppression  of  secondary  electrons 

at  the  collector  did  not  eliminate  or  even  significantly  reduce 

the  oscillation  of  positive  ions.  It  seems  therefore  that  secondary 

electrons  are  not  necessary  for  sustaining  the  oscillations  and 

that  some  other  feed-back  mechanism  must  also  exist.  It  was 
12 

suggested  by  Agdur  that  phenomena  in  the  gun  region  may  be 
significantly  related  to  the  oscillation  mechanism.  This  seems 
unlikely,  however,  since  the  transverse  oscillations  observed 
by  Mihran,  Agdur,  and  the  author  arc  identical,  even  though 
the  electron  guns  used  were  very  different  (Mihran's  gun  was 
gridded,  contained  a  planar  cathode,  and  was  shielded  from  the 
magnetic  field.  Agdur's  gun  also  was  gridded  and  contained  a 
planar  cathode,  hut  it,  was  completely  immersed  in  the  magnetic 
focusing  field.  The.  author's  gun  was  a  convergent-flow  Pierce 
gun  which  was  shielded  from  the  magnetic  field).  It  is  probable 
that,  the  transverse  oscillations  resulted  from  the  influence  of 
ions  or  electrons  (secondary  electrons  from  the  collector  or  from 
ions)  traveling  towards  the  electron  gun  combined  with  a  disturb¬ 
ance  of  the  profile  of  the  radial  potential  in  the  elecLron  beam 
similar  to  that  shown  in  Figure  L2. 

The  build-up  of  transverse  oscillations  might  be  accounted 
for  as  follows  Assume  first  that  the  electron  beam  was  not 
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oscillating  and  that  it  was  traveling  along  the  axis  of  the  magnetic- 
focusing  field.  If  the  ion  core  were  also  centered  on  the  axis  of 
the  solenoid,  then  nothing  would  happen;  however,  if  the  ion  core 
were  displaced  transversely,  as  is  shown  in  Figure  13a,  the 
electron  beam  would  follow  a  helical  path,  as  shown  in  Figure  13b. 
The  ion  core  would  not  remain  stationary,  but  would  move  toward 
the  center  of  the  electron  beam;  however,  when  the  ion  core  moved, 
the  electron  beam  would  assume  a  new  position  as  shown  in 
Figure  13c.  As  this  process  continued,  the  oscillations  of  the 
ion  core  would  increase  in  amplitude  because  the  electron  beam 
would  always  be  on  the  opposite  side  of  the  axis  of  the  solenoid 
from  the  ion  core  (see  the  appendix).  Now,  if  an  axial  potential 
of  the  appropriate  magnitude  increasing  in  the  plus  ■/.  direction 
were  applied  to  the  oscillating  ion  core,  it  would  move  in  the 
minus  •/,  direction  and  in  turn  cause  the  diameti  r  of  the  helical 
path  of  the  electron  beam  to  increase,  This,  of  course,  would 
further  enhance  the  ion  oscillations, etc .  The  oscillations  would 
stop  growing  when,  as  shown  in  Figure  13c,  the  ion  core  moved 
outside  the  electron  beam  so  that  ions  were  accelerated  towards 
the  walls  of  the  vacuum  chamber. 


The  axial  potential  resulted  from  the  penetration  of  the  cathodi 
potential  into  the  drift  region. 
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(C) 


Ion  core  tends  to  follow 
electron  beam  in  (b) 


In  many  cases  the  axial  potential  was  other  than  a. 
value  that  caused  continuous  enhancement  of  the  ion  oscilla¬ 
tions  so  that  beating  phenomena  occurred,  which  produced  the 
variety  of  amplitude  modulations  mentioned  earlier  in 
connection  with  the  transverse  oscillations. 

When  gridded  electron  guns  similar  to  those  used 
by  Agdur  and  Mihran  are  employed  to  generate  electron  beams 
for  plasma  devices,  ions  are  normally  not  accelerated  towards 
the  gun.  Electrons  freed  from  molecules  during  ionization, 
however,  are  accelerated  toward  the  gun,  and  these  can  cause 
a  feedback  similar  to  that  illustrated  in  rigure  13, 

C.  RADIAL  OSCILLATIONS 

Figure  14  gives  a  plot,  of  current  density  as  a  function 

of  time,  which  .shows  that  there  were  at  least  two  distinct  modes 

of  ion-induced  oscillations  in  the  electron  beam  at  a  beam  voltage 

-4 

of  5  kv,  a  magnetic  field  of  300  gauss,  and  a  pressure  of  3  x  10 
Tori'.  The  oscillations  shown  between  the  times  of  l|is  and  3.  Sps 
are  of  the  transverse  variety  described  in  the  previous  section. 
Tin;  fluctuations  prior  to  l|j.s  and  after  3.  5|is  result  from  oscilla¬ 
tion  of  the  beam  radius,  as  shown  in  Figure  15.  The  effect  of 
these  radial  oscillations  on  the  beam  voltage  and  on  the  magnetic 
field  was  found  to  be  veby  similar  to  that  of  the  transverse  oscilla  ¬ 
tions,  except  that  at  high  magnetic  fields  (  '>  500  gauss)  they 
di  sappeared.  The  frequency  was  found  to  be  independent  of  time 


and  press  ure  and  was  higher  than  the  frequency  of  the  trans¬ 
verse  osc  illations  by  a  factor  ranging  from  1,  4  at  the  beginning 
of  the  traxisverse  oscillations  to  2.  at  the  end,  as  can  be  seen 
in  Figure  14. 

B-y  comparing  the  predicted  frequency  of  the  transverse 
oscillatio  n 


with  that  of  the.  radial  oscillation,  Wj  r^,  it  is  seen  that  for  a 
small  ion  core,  (  c  <  <  d),  the  radial  frequency  should  be  V  2 
times  the  transverse  frequency.  This  compares  favorably 
with  the  value  of  1,  4.  When  c/d  is  not  small,  but  is  instead 
closer  tu  3/4,  it  is  apparent  that  the  radial  frequency  should 
be  about  l  times  the  transverse  frequency,  which  agrees  with 
the  uppo  r  limit  of  the  ratio  given. 

I  t  is  interesting  that,  since  the  frequency  of  the  radial 
oscillati  ons  was  independent  of  time  and  pressure,  the  plasma 
frequenc-.y  of  the  ions  Wj  was  also  conaLaul.  This  means  that 
changes  in  the  frequency  of  the  transverse  oscillations  as  a 
function  of  time  (Figure  H)  probably  resulted  from  an  increase 
in  the  raulius  of  the  ion  core  with  time. 

Sine  e  neither  Mihran  nor  Agilur  (nor  other  workers)  havi 
reported  radial  oscillations  of  the  type  shown  in  Figures  14  and 
it  is  probable  that  they  are  the  result  of  the  particular  gun 
geometary  used  for  the  nicasurm unis  made  in  this  study. 


-  lb 


It  is  shown  in  the  appendix  that  if  the  diameter  of  the  electron 
beam  oscillates  at  or  near  the  natural  frequency  of  the  radial 
oscillations  of  the  ion  core,  then  the  ion  oscillations  will  build 
up  exponentially  with  time.  In  examining  the  convergent  Pierce 
gun  for  possible  v/ays  in  which  oscillations  of  the  ion  core  could 
cause  the  diameter  of  the  electron  beam  to  oscillate,  it  was 
noted  that  a  density  modulation  of  the  ion  core  could  periodically 
alter  the  focusing  of  the  electron  beam  in  the  region  where  the 
beam  approached  its  minimum  diameter.  Indeed  a  density 
modulation  of  the  ion  core  is  to  be  expected  in  the  region  near 
the  gun  because  both  changes  in  the  radius  of  the  ion  core  and 
in  the  radius  of  the  electron  beam  affect  the  potential  along  the 
axis  of  the  beam,  which  in  turn  affects  the  velocities  of  the  ions 
and  their  density. 

The  absence  of  radial  oscillations  at  high,  levels  of  mag¬ 
netic  field  is  explained  by  there  being  enough  magnetic  flux 
linking  the  cathode  and  passing  through  the  region  where  the 
beam  approached  its  minimum  diameter  to  confine  the  beam. 

As  a  result,  no  feed-back  mechanism  existed,  and  the  oscil¬ 
lations  failed  to  grow,  The  initial  occurrence  of  radial  osc  il¬ 
lations  before  the  transverse  oscillations  built  up  were  the  result 
of  the  rapid  build-up  of  ions  in  the  beam  in  the  10'^  Tori'  range. 
Since  the  ion  build-up  time  was  comparable  to  the  period  of  the 
radial  oscillations,  the  ion  density  increased  to  a  large  value  be¬ 
fore  sufficient  time  elapsed  for  the  ions  to  move  radially  inward 
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into  the  potential  well  of  the  electron  beam.  The  ion  core 
then  decreased  in  diameter  and  oscillations  followed.  The 
reasons  for  the  oscillations  changing  from  radial  to  transverse 
and  back  to  radial  during  the  cathode  pulse  are  not  known. 

It  appears  that  radial  oscillations  of  the  electron  beam 
can  be  eliminated  by  using  magnetic  fields  well  above  the 
Brillouin  value  and  by  ensuring  that  the  magnetic  field  links 
the  cathode.  Apparently  there  is  no  convenient  way  to  prevent 
the  transverse  oscillations  of  the  beam  when  it  is  necessary 
to  operate  a  device  under  conditions  (of  pressure  and  time 
duration),  where  transverse  oscillations  normally  occur. 

The  only  possible  solution  seems  to  be  to  operate  at  a  pressure 
where  radial  oscillations  normally  occur  and  then  to  use  a 
magnetic  field  that  is  large  enough  to  prevent  the  transverse 
oscillations.  Tills  means  that  either  very  short  or  very  long 
pulses  should  he  used. 
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IV.  DISCUSSION  AND  CONCLUSIONS 


From  the  data  on  transverse  oscillations  and  from  the 
data  published  by  Mihran  and  Agdur,  it  can  be  concluded  that 
the  transverse  oscillations  increase  in  amplitude  with  time 
because  of  an  interaction  after  the  electrons  enter  the  magnetic 
focusing  field.  It  was  found  that, because  the  ion  core  disturbed 
the  potential  well  caused  by  the  electrons,  the  electron  beam 
followed  a  helical  path  from  the  gun  to  the  collector.  Ions  or 
electrons  traveling  towards  the  gun  altered  the  helical  path 
and  caused  the  oscillations  to  grow. 

Tlie  radial  oscillations  occurred  at  a  frequency  that  was 
1.4-2  times  the  transverse  frequency.  They  increased  in 
amplitude  with  time  because,  accompanying  the  oscillation  of 
the  ion  core,  there  was  a  density  modulation  of  the  ions  in  the 
core  that  affected  the  focusing  of  the  electron  beam  in  the  convcr 
gent  Pierce  gun  in  a  periodic  manner.  This  caused  the  diameter 
of  the  electron  beam  to  oscillate,  which  in  turn  caused  the  radial 
oscillations  of  the  ions  to  grow. 

It  is  concluded  that  the  radial  oscillations  can  be  eliminal 
by  using  a  strong  magnetic  field  that  links  the  cathode  and  beam- 
formation  regions.  It  appears,  however,  that  if  the  transverse 
oscillations  are  avoided  by  operating  at  high  pressure,  strong 
electron  cyclotron  oscillations  occur  which  may  counteract  any 
advantage  gained. 


V.  APPENDIX 


It  is  helpful  in  understanding  the  ion-induced  beam 
oscillations  if  the  work  of  Peterson  and  Puthoff^  is  extended 
to  include  the  effects  of  transverse  and  radial  oscillations  of 
the  electron  beam  at  or  near  the  plasma  frequency  of  the  ions. 

Shown  in  Figure  16  is  a  sketch  of  the  geometry  for  trans¬ 
verse  oscillation  to  be  analyzed.  The  potential  distribution  is 
found,  as  was  done  by  Peterson  and  Puthoff,  by  superposing  the 
potential  resulting  from  the  negative  and  positive  space  charges 
81  and  G-,  ,  and  a  space-charge -free  potential  0^  chosen  to 
make  the  sum  of  all  three  potentials  fil  the  boundary  condition  at 
the  anode  (whore  the  potential  equals  xcro). 


FIGURE  16.  Geometry  of  Transverse  Oscillations. 
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The  potential  is  given  by 


where 


2  .  _  "e 

V  “  ~ 


R  =  r2  +  -  2r£,  cos0  , 


and  is  the  electron  charge  density.  Then 


h  "  TT" 

o 


^ —  (r^-2  r  £,  cos  0) 


inside  the  negative  space-charge  region,  and  it  has  been  assumed 
that.  £,  «  r  .  Outside  the  negative  space-charge  region, 

v2  <|>j  =  0  ;  (A 

the rcforc 


,  p  /i  _  .. 

^1  ?  T—  ^  1  2  1,1  TT  “  2 -j  t:os°) 


In  a  similar  manner,  the  potential  in  the  positive  space -charge 
region  resulting  from  the  positive  space  charge  may  be  found 


— —  (r2  -  2r  6  cos  0) 
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where  it  lets  been  assumed  that  the  charge  density  of  the  ions 
equals  the  charge  density  of  the  electrons.  Outside  the  positive 
space -charge  region. 


'  PeC  r  6 

t)>^  =  -  ^  ^  (1  +  2  In  -  -  2  —  cost') 

o 

'  i  i 

(fit  Hie  anode  <|t(r,  0)  =  0  =  j  +  4*2  *  93  and  r  =  d  .  Therefore 
inside  the  anode 


p  a 

<t>3  =  -  g —  (1  +  2  In  ^  -  2  -2-j  r  cos 
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2  ~  r  eosO) 
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This  means  that  tin1  total  potential  is 


'I’l  '  (|>2  '  ‘I’d  ;  T 


k  \  ('  -  ?) '  '  (‘ 


l>  /  2  2 

>•  I  a  -  e 


2  C  2 


(A.  ID) 


"•'(•lore  the  force  equation  may  be  written  as 


<1Z  ,  .  (.  «  M  .  “l/f  /,  aZ\  , 
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Now  if  it  is  ossumod  Unit  f,  =  £  sinul  t  the  solution  is  of  Lho  fun 
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and 


(A.  121,) 

Thus  it  is  seen  that,  depending  on  the  value  of  to  ,  the  transverse 
oscillations  of  the  ions  may  be  amplitude  modulated  because  of 
sinwt  beating  with 


and  there  is  the  possibility  that  this  could  lead  to  amplitude  mod¬ 
ulation  of  the  transverse  oscillations  of  the  electron  beam. 

Shown  in  Figure  1  7  is  a  sketch  of  the  geometry  for  the 
radial  oscillations  to  be  analyzed,  Using  the  notation  of  Peterson 
and  Puthoff,  we  take  the  d-c  radius  of  the  plasma  column  as  being 

r  and  the  instantaneous  radius  as  r.  .  The  net  charge  per  unit 
o  i  ■ 

z 

length  inside  the  plasma  rore  rau^pfl  by  ions  is  rr  r  p.  ;  that 
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FIGURK  17.  Geometry  of  Radial  Oscillation  System. 

caused  by  olet  i  rons  is  -iri\  «  .  If  it  is  assumed  that  n  varies  as 

'  1  1  e  1  e 

p  “  P  I  p  CO  .4  tot  ,  (A 

e  1  eo  1 

win;  re  to  is  near  but  not  necessarily  equal  to  the  oscillation  frequency 
of  the  ions  in  the  plasma,  then  by  Gauss'  law  the  electric  field  at  the 
eel  in*  of  the  plasma  column  is 


i  O 

r  )  -  -= -  (p  t  +  P  cos  cot  -  *5 - — 

i'  2  c  '"eo  le~  2  r.  e 

n  1  O 
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therefore  the  force  equation  can  be  written  as 


”Z  R  +  fV"  ^eo  R  -  IT  +  R  Po-  COSwt) 

dt  o 
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where  r./r  =  R  .  When  r./r  is  near  unity,  1/R  =  2  -  R  ,  so 

l  o  l  o 

the  force  equation  can  be  rewritten  as 
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VALUES  OF  g 


FIGURF,  1R.  Iho-h  Stability  Chart  for 
Mathieu  Functions  of  Fractional  Order 


R  =  A  +  Be*12  sin  2a  +  Ce  sin  2e  t  D  sin  4z  ,  (A.  21) 

where  no  attempt  has  been  made  to  retain  the  phase  relationships 
between  the  sinusoidal  terms.  The  most  important  information 
contained  in  Equation  (A.  21)  is  that  as  long  as  a  remains  in  the 
unstable  region,  the  frequency  of  the  ion  oscillations  is  the  same 
as  the  frequency  at  which  the  electron  beam  oscillates,  and  that 
the  oscillations  build  up  exponentially  with  time.  This  means, 
of  course,  that  amplitude  modulation  similar  to  that  occurring 
with  transverse  oscillations  because  of  the  beating  of  the 
oscillation  frequency  of  the  electron  beam  with  the  natural 
frequency  of  oscillation  of  the  ion  core  should  not  occur  in  the 
case  of  radial  oscillations. 
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